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Evolutionary history of the sequestrate genus Rossbeevera 
(Boletaceae) reveals a new genus Turmalinea and highlights 
the utility of ITS minisatellite-like insertions for molecular 
identification 


T. Orihara’, T. Lebel?, Z.-W. Ge’, M.E. Smith’, N. Maekawa 


Abstract The sequestrate (truffle-like) basidiomycete genera Rossbeevera, Chamonixia, and Octaviania are closely 
related to the epigeous mushroom genera Leccinum and Leccinellum. In order to elucidate the properties and place- 
ment of several undescribed sequestrate taxa in the group and to reveal the evolutionary history of Rossbeevera and 
its allies, we conducted phylogenetic analyses based on three nuclear (ITS, nLSU, EF-1a) and two mitochondrial 
DNA loci (ATP6 and mtSSU) as well as precise morphological observations. Phylogenetic analyses of three nuclear 
loci suggest a complex evolutionary history with sequestrate fruiting bodies present in several clades, including a 
previously unrecognized sister clade to Rossbeevera. Here we propose a new sequestrate genus, Turmalinea, with 
four new species and one new subspecies as well as two new species of Rossbeevera. The three-locus nuclear 
phylogeny resolves species-level divergence within the Rossbeevera-Turmalinea lineage, whereas a separate 
phylogeny based on two mitochondrial genes corresponds to geographic distance within each species-level line- 
age and suggests incomplete lineage sorting (ILS) and gene introgression within several intraspecific lineages of 
Rossbeevera. Furthermore, topological incongruence among the three nuclear single-locus phylogenies suggests 
that ancient speciation within Rossbeevera probably involved considerable ILS. We also found an unusually long, 
minisatellite-like insertion within the ITS2 in all Rossbeevera and Turmalinea species. A barcode gap analysis 
demonstrates that the insertion is more informative for discrimination at various taxonomic levels than the rest of 
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the ITS region and could therefore serve as a unique molecular barcode for these genera. 
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INTRODUCTION 


In recent years, molecular systematics has become increasingly 
popular and useful for examining the evolution of fruitbody mor- 
phology within the fleshy fungi in Ascomycota and Basidiomy- 
cota. Most types of mushroom fruitbodies, including agaricoid, 
gasteroid, resupinate, cup-shaped, and coralloid forms, have 
arisen multiple times due to convergent evolution (e.g., Hibbett 
et al. 1997, Bruns et al. 1998, Percudani et al. 1999, Pine et al. 
1999, Peintner et al. 2001, Larsson et al. 2004, Hibbett 2004, 
Dentinger & McLaughlin 2006, Wilson et al. 2011, Miettinen et 
al. 2012). The sequestrate fruiting body form, which includes 
truffles, truffle-like fungi, and secotioid fungi, is characterized 
by fully or partially enclosed hymenia that render forcible spore 
discharge ineffective (Kendrick 1992). Sequestrate fungi have 
evolved dozens of times from other types of fruitbodies (e.g., 
Bruns et al. 1998, Percudani et al. 1999, Peintner et al. 2001, 
Læssøe & Hansen 2007, Lebel & Syme 2012; Lebel et al. 
2015). Previous statistical studies suggest that evolution of 
the sequestrate-gasteroid form from other fruitbody forms is 
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irreversible (Hibbett 2004, Wilson et al. 2011), and that this 
evolution may offer opportunities to adopt new ecological niches 
(Wilson et al. 2011). Many sequestrate fungi fruit belowground 
(hypogeous) and produce fruitbodies that are eaten by mam- 
mals, birds, or invertebrates that facilitate spore dispersal (Fo- 
gel & Trappe 1978, Maser et al. 1978, Claridge & May 1994, 
Danks 2011, Beever & Lebel 2014). Animal mycophagy may 
also help to disperse spores directly into similar habitats with 
compatible mycorrhizal host plants (Frank et al. 2009, Trappe et 
al. 2009). Since most sequestrate fungi form mycorrhizas with 
plants, this animal interaction should be beneficial for efficient 
spore dispersal. 


The Boletaceae is a family that primarily consists of soft-bodied, 
pileate-stipitate mushrooms with a poroid hymenophore (i.e., 
‘poletes’). In recent years, many new genera of boletes have 
been proposed based on multigene phylogenetic analyses of epi- 
geous members of the Boletaceae (e.g., Halling et al. 2012a, b, 
Arora & Frank 2014, Gelardi et al. 2014, Li et al. 2014, Wu et 
al. 2014, Zeng et al. 2014). However, the family includes many 
sequestrate lineages and the number of sequestrate genera 
also continues to increase as more taxa are included in mo- 
lecular phylogenies (Desjardin et al. 2008, 2009, Orihara et al. 
2010, Trappe et al. 2013). The sequestrate genus Rossbeevera 
was erected to accommodate Australasian and Asian species 
of Chamonixia (Lebel et al. 2012a, b) that are phylogenetically 
close to, but distinct from, Chamonixia s.str. Species of Ross- 
beevera have fruitbodies with a whitish peridium that discolours 
bluish green to indigo and ellipsoid to fusiform basidiospores 
with 3—5 longitudinal ridges. Phylogenetically, this genus is part 
of amonophyletic lineage that also includes the epigeous bolete 
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genera Leccinum and Leccinellum and the sequestrate genera 
Chamonixia and Octaviania (Lebel et al. 2012a, Orihara et al. 
2012b) (i.e., the ‘leccinoid clade’). Leccinum and Leccinellum 
are well known because many species in these genera are 
prized edibles collected throughout the Northern Hemisphere 
(but some Leccinum species can cause gastrointestinal dis- 
tress; Bessette et al. 2000). The leccinoid clade appears particu- 
larly rich in sequestrate fungi (> 40 spp.), even when compared 
with other groups of Boletales with many sequestrate taxa. A 
total of eight Rossbeevera species have thus far been reported 
from Australia (2 spp.), New Zealand (1 sp.), Southeast Asia 
(1 sp.), China (2 spp.), and Japan (2 spp.) (Lebel et al. 2012a, 
b, Orihara et al. 2012a). 


For several decades, a sequestrate fungus that is macroscopi- 
cally similar to Rossbeevera species but with pinkish fruitbodies 
and ellipsoid basidiospores with 6—10 longitudinal ridges has 
often been collected from western Honshu, Japan. Similarly, 
Rossbeevera-like fungi with different basidiospore morphology 
from Rossbeevera have also been formally and informally re- 
ported from other areas of Japan under alternative names (i.e., 
Gautieria sp. in Aoki (1978); Hymenogaster sp. in Yoshimi & 
Doi (1989)). All of these taxa have longitudinally ridged basidio- 
spores and in some cases exhibit a bluish discolouration that 
is similar to Rossbeevera species. In order to better assess 
the diversity of sequestrate fungi in the leccinoid clade, it is 
important to clarify the morphology and phylogeny of these 
poorly documented taxa. 


The aim of this paper is to provide well-resolved phylogenies 
based on three nuclear (ITS, nLSU, EF-1a) and two mitochon- 
drial (ATP6 and mtSSU) DNA loci (472 nucleotide sequences in 
total) and topological comparisons of the resultant phylogenies 
as well as precise morphological observations to enhance 
our understanding of evolution within the leccinoid clade. We 
include known species of Leccinum, Leccinellum, Chamonixia, 
Octaviania, and Rossbeevera as well as species that have not 
been formally studied in the past. Based on a combination of 
morphological and molecular characterization, we describe 
two new species of Rossbeevera as well as four new species 
and one new subspecies within the new genus Turmalinea. We 
also examine the utility of a highly divergent minisatellite-like 
insertion found within the ITS2 rDNA of Rossbeevera and Tur- 
malinea spp. for molecular delimitation within the Rossbeevera- 
Turmalinea lineage. 


MATERIALS AND METHODS 


Materials examined and morphological observations 


Fresh fruitbodies were collected throughout Japan and from 
Australia, New Zealand, and Yunnan Province, China. After 
morphological observation and DNA extraction, the fruitbodies 
were air-dried or freeze-dried for later examination. These 
specimens are deposited in Kanagawa Prefectural Museum 
of Natural History, Japan (KPM), National Museum of Nature 
and Science, Japan (TNS) and the Herbarium of Cryptogams, 
Kunming Institute of Botany, Chinese Academy of Sciences 
(KUN, with HKAS accession numbers). Herbarium specimens 
were obtained from KPM, TNS, the Oregon State University 
Herbarium, Oregon, USA (OSC), the Mycological Herbarium 
of Universita degli Studi di LAquila, L'Aquila, Italy (AQUI), the 
Western Australian Herbarium (PERTH), and the National 
Herbarium of Victoria, Royal Botanic Gardens Victoria, Aus- 
tralia (MEL). 


For microscopy hand-cut sections of fresh or dried specimens 
were mounted in water, 3 % KOH, lacto-glycerol, or 1 % phlox- 
ine B aqueous solution. To determine the amyloid reaction, 
dried material was stained with Melzer’s reagent. Basidiospore 


dimensions (e.g., range of spore length x spore width, length 
of hilar appendages) and their standard deviations (SD) were 
determined based on 50 measurements. The 95 % prediction 
intervals of basidiospore diameter are shown without paren- 
theses in taxonomic descriptions. Both endpoints of the spore 
dimensions are shown in parentheses, but when the value is the 
same as the 95 % prediction interval, only the latter is shown. 
Two additional spore features are shown; the length to width 
ratio (Q) and the hilar appendage to spore length ratio (HA/S; 
Orihara et al. 2012a). Measurements include the hilar append- 
age but not spore ornamentation or the pedicel. Basidium sizes 
are presented as the range of the lengths x the range of the 
widths. Scanning electron microscopy (SEM) was performed 
according to Yanaga et al. (2012). 


DNA extraction, PCR amplification and sequencing 


DNA was extracted from fresh or dried fruitbodies using the FTA 
Classic Card or Indicating FTA Cards (Whatman International 
Ltd, Maidstone, England) based on the manufacturer’s protocol 
for plants (www.whatman.com/References/WGI_1397_Plant- 
Poster_V6.pdf). Dried fruitbody tissues were soaked in 99.5 % 
ethanol for 2 min before DNA extraction according to Orihara 
et al. (2012a). PCR amplification of ITS, nLSU, and EF-1a fol- 
lowed Orihara et al. (2012b) and the protocol for mtSSU was the 
same as for nLSU. PCR amplification of ATP6 was performed 
according to Kretzer & Bruns (1999). PCR primers were ITS1F 
(Gardes & Bruns 1993) and ITS4 (White et al. 1990) or ITS4B 
(Gardes & Bruns 1993), ITS1F and ITS2, and/or ITS3 (White et 
al. 1990) and ITS4B for the ITS region; LROR and LR5 (Vilgalys 
& Hester 1990) for the nLSU; EF1-983F and either EF1-2218R 
or EF1-1567R (Rehner & Buckley 2005) for EF-1a; ATP6-3 
and ATP6-2 (Kretzer & Bruns 1999) for ATP6; and MS1 and 
MS2 (White et al. 1990) for mtSSU. Unidirectional sequencing 
of the PCR products in forward and reverse directions were 
completed according to Orihara et al. (2012b). Sequences were 
edited and assembled with Sequence Scanner v. 1.0 (Applied 
Biosystems, Foster City, California, USA), BioEdit v. 7.0.9 (Hall 
1999), SeaView v. 4 (Galtier et al. 1996) and Clustal X v. 1.83 
(Thompson et al. 1997). A total of 362 sequences were obtained 
and were deposited in the International Nucleotide Sequence 
Databases (INSD). 


Multigene phylogenetic analyses 


Specimens and sequences used for phylogenetic analyses are 
listed in Table 1. A total of 111 collections and 467 sequences 
(composed of 362 new sequences for this study and 105 se- 
quences from INSD) were included in the multigene datasets 
(nuclear dataset: ITS, nLSU, EF-1a; mitochondrial dataset: 
ATP6 and mtSSU). Identification of epigeous bolete taxa in the 
leccinoid clade that were included in the phylogenetic analyses 
were based not only on their sequence similarity to sequences 
in INSD but also on morphological descriptions by Den Bakker 
& Noordeloos (2005). Previous phylogenetic studies on boletoid 
mushrooms have shown that members of the leccinoid clade 
were closely related to Austroboletus spp., Tylopilus virens, 
Harrya chromapes (Binder & Hibbett 2006), Retiboletus griseus 
(Dentinger et al. 2010) or Spongiforma thailandica (Nuhn et 
al. 2013). Spongiforma thailandica is also closely related to 
Borofutus dhakanus and Porphyrellus species (Hosen et al. 
2012). We therefore included taxa from all of these groups 
in the analyses to see which is supported as a sister to the 
leccinoid clade based on different loci. Chalciporus piperatus 
and Buchwaldoboletus lignicola were strongly supported as the 
earliest diverging lineage in the Boletaceae (Binder & Hibbett 
2006, Dentinger et al. 2010, Halling et al. 2012a, b) and were 
therefore selected as outgroup taxa. Sequence alignments were 
performed using Clustal X with the data manually adjusted in 
SeaView v. 4. 
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KC552144 KC552218 


JN378406 
JN378436 
JN378437 
JN378432 
JN378440 
DQ029199 


JN378462 
JN378495 
JN378496 
JN378491 


JN257988 


KPM-NC 17751 


MEL 2341996 
MEL 2128484 


Japan, Hyogo Pref., Shiogahara 


Octaviania nonae 


KC552156 KC552231 
KC552157 + KC552232 


KC552004 


KC552005 
JQ619180 


KPM-NC 17819 


MES270 


KC552236 


KC552161 


JN378498 


DQ200917 AY700189 


HKAS 41551 (AFTOL-ID450) 


CD567 


KC552102 


KC552061 
KC552016 KC552058 KC552099 


MEL 2053830 


KC552017 KC552059 KC552100 


MEL 2300518 
MEL 2265009 


KC552015 KC552057 KC552098 


KC552014 KC552056 KC552097 


MEL 2233764 
MEL 2305143 


MB 03-118 


Australia, Tasmania, Mount Field, Mt Field Nat. Park. 


Octaviania tasmanica 


Australia, NSW, Southern Tablelands, off Nungatta Rd. 


Japan, Okinawa Pref., Ishigaki Isl., Mt Omoto 
USA, Maine, Tunk Lake, off route 182 


China 


Octaviania yaeyamaensis 


Octaviania zelleri 


Aureoboletus thibetanus 


Australia, Victoria, Aire River Crossing, Otways National Park 


Austroboletus occidentalis 


Austroboletus sp. 1 


Australia, Victoria, On Quarry Track off Wannon Diversion Track about midway between 


Wannon Diversion Track and quarry 


Australia, Tasmania, Peter Murrell Nature Reserve. Site HEB. (55G, 524413, 5238290) 


Australia, Victoria, Colac Otway Shire, Carlisle State Park, Cricket Pitch Track, 


Austroboletus sp. 1 


Austroboletus sp. 2 


7.5 km W Gellibrand 
Australia, Western Australia, Walpole - Nornalup National Park, the Knoll lower walk 


New Caledonia, Southern Province, Thio to Nagety Road, 26 km north of Thio 


USA, MA 


Austroboletus sp. 3 


KC552018 KC552060 KC552101 


AY789082 


Austroboletus sp. 4 


AY879116 


AY684158 


Boletellus projectellus 
Borofutus dhakanus 


JQ928617 JQ928575 


JQ928607 


HKAS 73792 


Pul1 


Bangladesh, Dhaka Division, Gazipur, Bhawal National Park 


Germany, Maindreieck 


JQ327040 


JQ326997 


Buchwaldoboletus lignicola 
Chalciporus piperatus 
Harrya chromapes 


DQ534648 GU187690 


JN378517 


MB 04-001 


USA, MA, Rutland State Park 


KC552173 KC552249 


JN378457 


KC552019 


KPM-NC 17835 
MB 97-023 


Japan, Tottori Pref., Yazu-cho 
Germany, Walhalla, Bavaria 


DQ534563 DQ534643 GU187734 


Porphyrellus porphyrosporus 


Retiboletus griseus 


HQ161858 KF030414 


BD210 (LSU), Both sn (EF-1a) 


USA, VA, George Washington National Forest, Shenandoah Mountain (LSU sequence); 


NY (EF-1a sequence) 


China 


JQ928626 JQ928580 


JQ928613 


HKAS59460 


Retiboletus aff. griseus 


Rhodactina incarnata 


KC552020* 


OSC 130505 


Thailand, Chang Mai, Sanpatong District, Mae Wang Concervation Forest, Sanpatong- 


Ban Guard Road 


EU685108 KF030436 


EU685113 


DED 7873 


Thailand, Nakorn Nayok Province, Khao Yai National Park 


Japan, Kanagawa Pref., Odawara-shi, Iryuda 


Spongiforma thailandica 


Tylopilus virens 


KC552062 KC552103 KC552174 KC552250 


KPM-NC 18054 


For testing topological congruence among different loci, we 
directly compared statistically supported topologies (BS > 75 %) 
in ML trees constructed with RAxML v. 7.2.6 (Stamatakis 2006), 
setting the number of bootstrap replicates to 1 000. Introns of 
EF-1a were carefully adjusted in SeaView and also included in 
the analyses. Approximately 80 % of nucleotides were invaria- 
ble in the mtSSU dataset and we therefore included only 106 bp 
of variable sites in the phylogenetic analyses. 


In the ITS dataset, a large part of the ITS1 and ITS2 regions 
were too divergent to be aligned across the leccinoid clade and 
the allied taxa. Therefore, we first tested whether to include 
unambiguously aligned sites using the Gblocks algorithm (Cas- 
tresana 2000) in SeaView v. 4. This allowed smaller final 
blocks, gap positions within the final blocks, and less strict 
flanking positions. The resultant dataset was 263 bp in length 
and many phylogenetically informative positions that contained 
gaps were excluded from the analysis based on the results pro- 
duced by Gblocks. Consequently, the ML tree generated from 
the dataset by RAxML was poorly resolved with evidence of 
long-branch attraction of B. dhakanus and S. thailandica within 
the Leccinum clade. Because of these inconsistent results we 
abandoned the Gblocks dataset and manually and carefully 
selected phylogenetically informative sites within ITS using 
SeaView, allowing inclusion of the sites that contained gaps 
(up to approximately 40 % of the individuals at any given site). 
The final ITS dataset for the multi-locus analyses was 462 bp 
in length, and the resultant ML tree had much higher resolution 
without obvious evidence of long-branch attraction. 


Subsequently, two multigene datasets that consisted of three 
nuclear (i.e., ITS, nLSU, and EF-1a) and two mitochondrial 
(i.e., ATP6 and mtSSU) DNA loci, respectively, were created to 
examine the phylogeny of the leccinoid clade and to check for 
topological incompatibility between nuclear and mitochondrial 
phylogenies. In both datasets, we included only specimens with 
sequences for at least two loci (i.e., specimens represented by 
only ITS were omitted in the nuclear dataset; Table 1). Gaps 
were treated as ‘missing’ data for all analyses. The nuclear and 
mitochondrial combined datasets are deposited in TreeBASE 
(accession URL: http://purl.org/phylo/treebase/phylows/study/ 
TB2:S15097). 


For each dataset, Bayesian analysis was conducted with Mr- 
Bayes 3.2 (Ronquist & Huelsenbeck 2003). The combined 
dataset was partitioned by regions for the non-coding se- 
quences and by codons for the coding sequences, and the 
best-fit likelihood models were estimated for the resultant 10 
partitions with MrModeltest 2.3 (Nylander 2004). The GTR+I+G 
model was used for all the nuclear partitions (i.e., 1TS1-5.8S- 
ITS2, nLSU, and the 1st, 2nd, and 3rd codons and introns of 
EF-1a) and the 1st and 3rd codons of ATP6, the GTR+G model 
for mtSSU, and the F81 model for the 2nd codon of ATP6. 
Bayesian posterior probabilities (PP) were approximated by the 
Metropolis-coupled Markov chain Monte Carlo method (Geyer 
1991). Two parallel runs were conducted with one cold and 
seven heated chains each for 10M generations, starting with 
a random tree. The parameter for temperature of the seven 
heated chains in both runs was set to 0.10. Trees were saved 
to a file every 1 000th generation. We determined that the 
two runs reached convergence when the average SD of split 
frequencies (ASDSF) continuously dropped below 0.01. We 
further verified the convergence by checking that the effective 
sample size (ESS) of each resulting statistic was sufficiently 
large (> 200) using Tracer v. 1.6 (Rambaut & Drummond 2009). 
Trees obtained before reaching convergence were discarded as 
the burn-in, and the remaining trees were used to calculate a 
50 % majority consensus topology and to determine PP values 
for individual branches. 
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The same combined datasets were also analysed by the maxi- 
mum likelihood (ML) method with raxmIGUI 1.3 (Silvestro & 
Michalak 2012), which includes executable files of RAxML 7.4.2 
(Stamatakis 2006). The best-fit ML tree was inferred under the 
GTR+I+G model. The datasets were partitioned in the same 
way as in the Bayesian analysis so that different a-shape pa- 
rameters, GTR rates, and empirical base frequencies could be 
assigned to each partition. To check statistical support for the 
resultant tree topology, the rapid bootstrap (BS) option was used 
under the automatically assigned GTR+CAT model, setting the 
number of replicates to 1 000. 


Species tree analyses 


Since several topological conflicts were found among individual 
single-locus nuclear phylogenies, we conducted species tree 
estimation from gene trees using *BEAST program (Heled & 
Drummond 2010) in the software package BEAST 2 v. 2.1.3 
(Bouckaert et al. 2014) to see whether the multi-locus combined 
phylogeny conflicts with the species phylogeny based on a 
multispecies coalescent model. *BEAST estimates a Bayes- 
ian species tree from the simultaneously generated, individual 
gene trees, and reconstructs a species phylogeny that reflects 
the effect of incomplete lineage sorting (ILS) or deep coales- 
cence (Heled & Drummond 2010). In the *BEAST analysis 
we excluded species with any missing loci from the analysis 
because they could negatively effect the MCMC convergence 
(Castillo-Ramirez et al. 2010). Thus, we were not able to in- 
clude Austroboletus cf. novae-zelandiae, Buchwaldoboletus 
lignicola, Chalciporus piperatus, Leccinum ‘albellum’, Leccinum 
sp., L. versipelle, Octaviania arbucalensis, O. asterosperma, 
O. zelleri, Retiboletus griseus, ‘Turmalinea chrysocarpa’, and 
Tylopilus virens. For the other species we included specimens 
with sequences of all the three loci. Consequently, the con- 
catenated nuclear datasets consisted of 78 sequences, which 
were grouped into 43 species according to the results of the 
ITS, nLSU, and EF-1a analyses. The dataset was partitioned 
into ITS, nLSU, and EF-1a regions. 


XML command files for executing *BEAST analysis was created 
using BEAUti in the BEAST 2 package. Substitution models for 
each locus were estimated with MrModeltest: the HKY+I+G 
model was assigned to the non-partitioned EF-1a dataset 
whereas the GTR+I+G model was used for the ITS and nLSU 
datasets. The number of rate categories for the discrete gamma 
model was set to four according to Yang (1994), and ‘piece- 
wise linear and constant root’ was applied to the species tree 
population size model. In the Bayesian MCMC prior settings, 
gamma distribution was applied to mutation rates and clock 
rates for each dataset. Both the analyses were run for 200M 
generations, storing an inferred tree in every 10 000 genera- 
tions. The other settings were left at defaults. The runs were 
monitored using Tracer v. 1.6 to verify the convergence of the 
MCMC by checking that the ESS of each resulting statistic was 
sufficiently large (> 100). The output data were summarized 
using TreeAnnotater in the BEAST 2 package. The resulting 
trees were viewed with FigTree v. 1.4. (http://tree.bio.ed.ac. 
uk/software/figtree/). 


Detection and characterization of insertion sequences 
within ITS2 


The full-length ITS dataset (3’ end of SSU-ITS1-5.8S-ITS2-5’ 
end of LSU) of the Rossbeevera and Turmalinea species 
consisted of 51 sequences. The dataset was not properly 
aligned by the automatic alignment with ClustalX v. 1.83 or 
with Q-INS-I option in MAFFT v. 7 (Katoh & Standley 2013) 
due to highly divergent sequences with indels and unusually 
long insertions within ITS. We therefore aligned the dataset 
manually with SeaView. The dataset is deposited in TreeBASE 


(accession URL: http://purl.org/phylo/treebase/phylows/study/ 
TB2:S15097). Internal, discontinuous repeats within the ITS 
region were detected with the aid of dot plot matrices generated 
with the BLASTN search using the default settings (Altschul et 
al. 1997). The search of more-or-less regular, continuous, ITS 
tandem repeats was implemented with Tandem Repeats Finder 
software v. 4.04 (Benson 1999) under the following parameter 
settings: match = 2, mismatch = 3, indel = 5, minimum align- 
ment score to report repeat = 50, and maximum period size = 
300. In one case the analysis detected both a single nucleotide 
motif and the repeated pattern that consisted of the same core 
motif and its successive duplication as ‘consensus patterns’ 
(i.e., both ‘TTTGCA and ‘TTTGCATTTGCA)). In that case, we 
counted only the single nucleotide motif as a detected consen- 
sus pattern (i.e., ‘TTTGCA in the above example). After careful 
alignment of the dataset, the boundaries of the insertions were 
determined based on the results of BLASTN searches and 
Tandem Repeats Finder analyses. 


To assess the usefulness of the ITS2 insertion for taxonomic 
identification, we implemented barcode gap analysis using the 
web version of the ABGD (Automatic Barcode Gap Discovery) 
software, which was primarily designed to delimit species from 
unexplored groups of metazoans (Puillandre et al. 2012). ABGD 
detects gaps between the infraspecific and interspecific diver- 
gence based on a matrix of pairwise distances and groups the 
input sequences into several hypothetical species. The larg- 
est gap is initially detected with various prior limits of maximal 
distance of infraspecific divergence (P), so that the dataset is 
partitioned according to the gap. The gap detection is iterated 
within each recursive partition until ‘large gaps’ are no longer 
found under the parameter P and X, where X is a proxy for 
minimal gap width (Puillandre et al. 2012). The analysis was per- 
formed using the default settings (i.e., Pu = 0.001, Pax = 0.1, 
X = 1.5) except the number of steps for P and for a distance- 
distribution histogram were set to 100 and 60, respectively. 
The web version of ABGD provides three distance models for 
the analysis: Jukes-Cantor, the Kimura’s two-parameter, and 
simple distances. Accordingly, we chose the Kimura’s two- 
parameter model (Kimura 1980; Ti/Tv = 0.9059 in the ITS2 in- 
sertion dataset; 1.2955 in the residual ITS dataset excluded 
the insertion sequences), which had a higher likelinood score 
than the Jukes-Cantor model based on model estimations in 
MrModeltest 2.3. In addition, ML analysis of the ITS2 inser- 
tion dataset was implemented using RAxML under the GTR+1 
model, which was chosen as the best model in MrModeltest 
2.3, and rapid bootstrap analysis was carried out in the same 
settings as those of combined datasets described above. 


RESULTS 


Topological incongruence among different single-locus 
datasets 


Prior to multi-locus analyses, we first examined topological 
congruency among three nuclear-locus trees (ITS, nLSU, and 
EF-1a) and among two mitochondrial-locus trees (ATP6 and 
mtSSU) separately. The resulting tree topology in the ITS, 
nLSU, and EF-1a datasets showed conflicts in the branching 
order within the Rossbeevera clade (Fig. 1; entire trees not 
shown). Maddison & Knowles (2006) suggest that including 
more individuals and more loci reduces the effect of incom- 
plete lineage sorting and contributes to reconstructing a more 
accurate species tree. Since all of the species in this clade 
contain multiple individuals, we concatenated the datasets for 
the following analyses and carefully compared the resultant 
multi-locus phylogeny with a species tree inferred from the 
three single-locus Bayesian phylogenies by using *BEAST. In 
contrast to the nuclear dataset, there was no obvious conflict 
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ITS rDNA 


LSU rDNA 


Australasian 
Rossbeevera spp. 


‘R. paracyanea’ 


R. griseovelutina 


R. eucyanea 
(Honsyu & Kyusyu) 


R. eucyanea 
(Okinawa & Amami-oshima Isls) 


| ‘R. cryptocyanea’ 


I R. yunnanensis 


0.05 substitutions / site 


‘R. paracyanea’ 


New lineage (‘Turmalinea’) 
I R. yunnanensis 


0.05 substitutions / site 


EF-1qa 


Australasian 
Rossbeevera spp. 


Australasian 
Rossbeevera spp. 


R. eucyanea 
R. griseovelutina 


| ‘R. cryptocyanea’ 
| ‘R. cryptocyanea’ 


R. griseovelutina 4 , 
g R. paracyanea 


R. eucyanea 


I R. yunnanensis 


0.05 substitutions / site 


Fig.1 Comparison of maximum likelihood (ML) tree topologies of three nuclear-locus phylogenies (ITS, nLSU, and EF-1a) of Rossbeevera. Branches sup- 
ported by RAxML bootstrap (BS) values > 75 % (1 000 replicates) are shown as thickened lines. Names with dashes represent the lineages that have not 


been reported in previous studies. 


between the mitochondrial ATP6 and mtSSU phylogenies, partly 
due to poor resolution of the mtSSU phylogeny. 


Although most of the species represented by multiple individuals 
formed monophyletic clades in the three nuclear datasets, we 
found topological discordance in the relationship in Octaviania 
subg. Parcaea, which resulted in collapse of the monophyly of 
the species-level clades. That is, in the nLSU phylogeny the 
O. nonae sequences showed paraphyly with O. celatifilia. This 
topology is documented in Orihara et al. (2012b) and has only a 
minor influence on the phylogenetic relationships of Rossbee- 
vera and its relatives, which are the focus of this study. Thus, 
we retained the sequences in the combined dataset. 


Species tree and multi-locus combined tree based on the 
nuclear datasets 


Since topological discordance among the three nuclear phylo- 
genies was observed in the Rossbeevera clade (Fig. 1), we first 
inferred a species tree from Bayesian ITS, nLSU, and EF-1a 
‘gene’ trees using *BEAST and subsequently compared the 
topology with the multi-gene tree based on the concatenated 
dataset. In the *BEAST analysis, all the parameters reached 
plateau before 5M generations. We therefore discarded the 
first 500 species trees as burn-in and summarized the remain- 
ing 19 501 trees for the PP calculation. ESS values of all the 
parameters became sufficiently large (> 100) after 200M gene- 
rations. The average Log Likelihood (In L) of the sampled spe- 
cies trees was 19 907.124. 


The combined dataset of three nuclear loci (i.e., ITS and LSU 
rDNA, and EF-1a) consisted of 110 individuals and 2 586 base 
pairs. In the Bayesian inference, the ASDSF of the two parallel 
MCMC runs dropped constantly below 0.01 after c. 2.15M ge- 
nerations. Thus, the first 2 150 trees in each run were discarded 
as burn-in according to the ASDSF value, and the remaining 
15 702 trees (sampled from 7.85M generations) were used 
to calculate a 50 % majority consensus tree and determine 
PPs. ESSs of all the parameters were sufficiently large. Total 
arithmetic mean and harmonic mean of Likelihoods (In L) were 
—23 076.11 and —23 171.84, respectively. The maximum likeli- 


hood analysis resulted in one ML tree (In L = —22 947.590387). 
Overall tree topologies between the Bayesian and ML trees 
were almost identical, and no topological conflict was observed. 
The resultant topology recovered a strongly supported lecci- 
noid clade (BS = 100 %, PP = 1.00) that included Chamonixia, 
Leccinum, Leccinellum, Octaviania and its three subgenera, 
Rossbeevera, and a novel sequestrate lineage sister to Ross- 
beevera (Fig. 3) (i.e., Turmalinea). 


Although there were some taxa omitted in the *BEAST species 
tree, overall topologies of the species tree and three-locus 
combined tree were similar (Fig. 2, 3). There was no strong 
topological conflict found between the two phylogenies. Both 
phylogenies showed that the novel Turmalinea clade was mono- 
phyletic and sister to Rossbeevera spp. In the *BEAST spe- 
cies tree the new lineage ‘R. cryptocyanea’ was shown to be 
a sister to R. eucyanea with strong statistical support despite 
the topological inconsistency among the single nuclear gene 
trees. However, the other infrageneric conflicting topologies 
among the gene trees were not resolved in the species tree 
analysis although they were resolved in the multigene combined 
analyses with relatively high statistical support. We therefore 
refrain from highlighting these infrageneric relationships here 
and will discuss it in more detail below. 


In the Bayesian and ML combined phylogenies, every genus 
and subgenus in the leccinoid clade except Leccinellum was 
strongly supported as monophyletic by both BS and PP (Fig. 3). 
The novel Turmalinea clade contains four undescribed se- 
questrate species and two subspecies from Japan and China. 
All of the lineages were strongly supported by both BS and 
PP (BS 2 89 %, PP = 1.00). In the Rossbeevera clade, two 
unknown species-level lineages were nested within the clade; 
Rossbeevera cryptocyanea sp. nov. and R. paracyanea sp. nov., 
both of which are morphologically similar to R. eucyanea and 
are distributed in Japan. 


The sequestrate genus Chamonixia was recovered as the earli- 
est diverging lineage within the leccinoid clade with moderate 
statistical supports (Fig. 3). The sister lineage to the leccinoid 
clade was not resolved in our analyses. A sequestrate species, 
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Aureoboletus thibetanus 
Retiboletus aff. griseus CHINA 
Harrya chromapes 
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Borofutus dhakanus 
Porphyrellus porphyrosporus 


Fig. 2 Bayesian species tree inferred based on three nuclear loci (ITS, nLSU, and EF-1a) using BEAST. A total of 78 individuals with sequence data of all 
three loci were grouped into 43 terminal species-level taxa. Bayesian posterior probabilities (PP) are indicated above or below branches. Branches strongly 


supported by PP are highlighted as thickened lines. 


Spongiforma thailandica, that was formerly supported to form a 
monophyletic clade with R. griseus and Leccinum spp. (Nuhn 
et al. 2013), formed a strongly supported monophyletic clade 
with Borofutus dhakanus and Porphyrellus porphyrosporus as 
shown in Hosen et al. (2012). Consequently, a recently pro- 
posed subfamily Leccinoideae, which includes Retiboletus, 
Borofutus, and Spongiforma as well as the genera in the lecci- 
noid clade (Wu et al. 2014), was not recovered in our nuclear 
multi-locus analyses. 


Mitochondrial multi-locus phylogeny of the Leccinoid clade 


The two-locus mitochondrial dataset (i.e., ATP6 and SSU 
mtDNA) consisted of 85 specimens and was 840 bp in length. 
In the Bayesian inference, the ASDAF of the two parallel MCMC 
runs dropped constantly below 0.01 after c. 1.77M generations. 
Thus, the first 1 770 trees in each run were discarded as burn-in, 
and the remaining 16 462 trees (sampled from 8.23M genera- 
tions) were used to calculate a 50 % majority consensus tree 
and determine PPs. ESS values of all the parameters were 
sufficiently large (> 200). Total arithmetic mean and harmonic 
mean of Likelihoods (In L) were —3 541.11 and —3 614.29, 
respectively. The In L of the resultant ML tree in the RAxML 


analysis was —3 431.324125. Monophyly of Leccinum, Ross- 
beevera, the Turmalinea clade, and the three Octaviania sub- 
genera were recovered with varying statistical support (Fig. 4). 
However, most basal topologies, including relationships of the 
three Octaviania subgenera, were not resolved or were poorly 
supported. 


Topological comparison between the nuclear and mito- 
chondrial phylogenies 


We examined topological differences between nuclear and mito- 
chondrial phylogenies (Fig. 3, 4) and compared the resultant 
topologies of the Rossbeevera-Turmalinea lineage to determine 
whether there were any discrepancies in patterns of divergence 
(Fig. 5). The nuclear phylogeny based on the combined data- 
set of ITS, LSU, and EF-1a had high resolution and generally 
showed infraspecific divergence with strong statistical support 
regardless of the geographic distribution of the samples. For 
example, sequences of R. griseovelutina, ‘R. paracyanea’ and 
‘R. cryptocyanea’ diverged into two strongly supported infra- 
specific lineages (lineages 1 and 2 in each species; Fig. 5). 
However, we did not find evidence of any obvious biogeographic 
disjunction within these species nor did we detect distinct mor- 
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Fig. 3 Bayesian 50 % majority-rule consensus tree of the nuclear three-locus dataset (ITS, nLSU, and EF-1a) of Rossbeevera and allied genera (i.e., the 
leccinoid clade). Bayesian PP and RAxML bootstrap (BS) values (1 000 replicates) are indicated above or below branches or at nodes as PP/BS. Values of 
PP < 0.90 or BS < 50 % are not shown. Branches supported by both PP 2 0.98 and BS 2 70 % are depicted as thickened black lines. Branches supported 
by both PP 2 0.98 and BS < 70 % are shown as thickened grey lines. Incompatible topologies between the combined phylogeny and nuclear single-locus 
phylogenies are indicated by red asterisks. Names of new taxa and new lineages are coloured in red. Chalciporus piperatus and Buchwaldobobletus lignicola 
were used as outgroups. 
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Fig. 4 Cladogram of Bayesian 50 % majority-rule consensus phylogeny of the mitochondrial two-locus dataset (ATP6 and mtSSU). Bayesian PP and RAxML 
BS are indicated at nodes as PP/BS. Values of PP < 0.90 or BS < 50 % are not shown. Branches supported by both PP > 0.98 and BS > 70 % are depicted 
as thickened black lines. Branches supported by either PP > 0.98 or BS > 70 % are shown as thickened grey lines. Names of new taxa and new lineages are 


coloured in red. 
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Fig.5 Comparison of nuclear three-locus (ITS, nLSU, and EF-1a; left) and mitochondrial two-locus (ATP6 and mtSSU; right) ML topologies of the Rossbeevera 
and Turmalinea clade. Bayesian PP and RAxML BS are indicated at nodes as PP/BS. Values of PP < 0.90 or BS < 50 % are not shown. The collection locali- 
ties are indicated on branches of either the nuclear or mitochondrial phylogenies. Location data are as follows: TAS = Tasmania, VIC = Victoria WA = Western 
Australia, NSW = New South Wales, NNZ = North Island, New Zealand, SNZ = South Island, New Zealand, W = Western, E = Eastern. 


phological differences between the infraspecific lineages. Little 
divergence was recognized within Australasian Rossbeevera 
lineages, at least at the species level, despite the different 
morphologies and distributions of the described species. On 
the other hand, the mitochondrial phylogeny had relatively 
poor resolution and also generally reflected differences in dis- 
tribution within each species shown in the nuclear phylogeny. 
Within the terminal Australasian Rossbeevera clade, we noted 
distinct divergence between the nuclear and the mitochondrial 
datasets. The mitochondrial topology of the clade was con- 
cordant with distributional patterns rather than morphological 
circumscription proposed by Lebel et al. (2012a) although some 
branches were not statistically supported (Fig. 5). For example, 
specimens that were morphologically identified as R. vittatispora 
were remarkably divergent in the mitochondrial dataset, but 
clearly reflected geographical proximity between Tasmania and 
Victoria. In the mitochondrial phylogeny, one individual of the 
‘R. paracyanea’ lineage 1 (KPM-NC 18087) was divergent from 
the ‘R. paracyanea’ lineage 2, whereas another individual of the 
‘R. paracyanea’ lineage 1 (i.e., KPM-NC 17847) was identical 
to the ‘R. paracyanea’ lineage 2. A very similar pattern was 
observed between ‘R. cryptocyanea’ lineages 1 and 2, although 
they also exhibited geographical proximity among the speci- 
mens in the mitochondrial phylogeny (i.e., Kakeroma Island is 


separated from Amami-oshima Island only by a narrow strait 
and they belong to the same group of islands, whereas Ishigaki 
Island is about 650 km distant from Amami-oshima). 


Characterization of a minisatellite-like insertion within ITS2 


We found that species of Rossbeevera and the Turmalinea 
clade have an unusually long insertion with indistinct tandem 
repeats within their ITS2 rDNA region. The inserted sequences 
ranged in size from 199 to 603 bp (Table 2). The Tandem Re- 
peats Finder analysis detected one or two continuous tandem 
repeats in the ITS2 of 12 of 16 lineages of the clade (Table 2). 
Most of the detected tandem repeats contained a common 
GAGTGAAAGTG motif or similar derivative sequences. This 
core pattern and the derivatives scattered within insertions of 
every lineage in the group. In addition, one terminal lineage 
of the Turmalinea clade (i.e., ‘T. persicina’) had another highly 
conserved core motif, AATTTAATTAT TGGTGTGGAAGCTT- 
GATTGTAAA, which is 34 bp in size and appeared three times 
just after the 3'- end of the GAGTGAAAGTG core pattern in 
the insertion. Turmalinea persicina is divided into two slightly 
divergent infraspecific groups (Fig. 6c). The copy number of 
the continuous motifs detected by the Tandem Repeats Finder 
analysis ranged from 1.9 to 5.9 (Table 2). 
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Fig. 6 Barcode gap analysis and ML phylogeny of minisatellite-like insertion within the nuclear ITS2 rDNA of Rossbeevera and the Turmalinea clade. — 
a. Results of the Automatic Barcode Gap Discovery (ABGD) analysis of the minisatellite-like insertion dataset. The dataset is composed of 51 sequences and 
is 838 nucleotides in length. Histogram of genetic distances between sequences and cumulative frequency of the distance value (ranked value; indicated as 
the blue line) in the distance matrix are shown. The ‘barcode gaps’ are represented as horizontal lines on the ranked value. — b. Results of the ABGD analysis 
using the ITS dataset excluding the minisatellite-like regions. The dataset is composed of the same individuals as ‘a’ and is 668 nucleotides in length. — c. ML 
phylogeny based on the same insertion dataset as ‘a’. 10 species-level groups recovered with prior infraspecific divergence 0.019 < P < 0.038 in ABGD are 
designated onto the phylogeny as [G1] — [G10]. These partitions correspond to the boundary of interspecific divergence designated in ‘a’. Two infraspecific 
groups of T. persicina are identical to those in Table 1. The RAxML bootstrap (BS) value (1 000 replicates) is indicated above or below branches (only BS 
values 2 50 % are shown). 


186 


Persoonia — Volume 37, 2016 


The minisatellite-like insertions found in every species of the 
Rossbeevera-Turmalinea lineage were highly conserved within 
each species (and infraspecific lineage) but highly divergent 
between species or genera (Table 2). Furthermore, the barcode 
gap analysis showed that the minisatellite sequence region can 
be used to delimit the boundary between intergeneric-level, 
interspecific-level, and infraspecific-level divergence more dis- 
tinctly than the rest of the ITS region (Fig. 6a, b). The boundary 
between infra- and interspecific-level divergence is supposed to 
be present at 0.0025 < D < 0.003, where Dis a genetic distance 
value of pairwise distances in a dataset (Fig. 6a). The minisatel- 
lite dataset was both primarily and recursively partitioned into 
the ten species-level clusters already represented in Fig. 3, 5 
with 0.019 < P < 0.038 in the recursive partitions, where P is a 
prior limit of maximal distance of infraspecific divergence (Fig. 
6c). As in the other nuclear analyses, Australasian Rossbeevera 
spp. were not very divergent and this species-complex was not 
partitioned at the species-level genetic boundary estimated by 
the barcode gap analysis. Also, there was a distinct difference 
in interspecific genetic distance within Rossbeevera and within 
Turmalinea in the minisatellite dataset, while the difference was 
unclear in the rest of the ITS region (Fig. 6a, b). 


The ITS2 minisatellite ML phylogeny mostly recovered mono- 
phyly of the Rossbeevera and Turmalinea species repre- 
sented in the other phylogenies in this study (Fig. 6c; In L = 
—2 615.994733). The monophyly of R. griseovelutina was not 
supported in the ML phylogeny of the minisatellite-like inser- 
tion within the nuclear ITS2 rDNA, probably due to the effect 
of the long genetic distance between the Rossbeevera and 
Turmalinea clades. The infraspecific divergence within ‘R. para- 
cyanea’, R. griseovelutina, ‘T. persicina’, and ‘T. mesomorpha’ 
was strongly supported in the analysis. Contrary to the pattern 
found in species of the Rossbeevera-Turmalinea lineage, the 
sister taxon of this group, Leccinellum griseum (= L. pseudosca- 
brum), had only a single 111 bp insertion in the ITS 2 region and 
Tandem Repeats Finder analysis did not detect any repeating 
motifs in this DNA fragment. 


Taxonomy 


The phylogenetic analyses conducted in this study consistently 
recovered a distinct yet previously unrecognized lineage that 
was sister to the monophyletic genus Rossbeevera. Morphologi- 
cally, members of the two clades are readily distinguishable by 
differences in the basidiospores and in the sterile base of the 
fruiting body. The unique, minisatellite-like insertions within ITS2 
nuc-rDNA were highly divergent between the two clades as well 
as between the different species-level lineages (Fig. 6). Based 
on both molecular and morphological data, we propose the new 
genus Turmalinea for this lineage. This group is comprised 
of three new Japanese species (T. persicina, T. yuwanensis, 
and T. mesomorpha), one subspecies (T. mesomorpha subsp. 
sordida) and one Chinese species (T. chrysocarpa). Within the 
genus Rossbeevera we propose two new species from Japan, 
R. paracyanea and R. cryptocyanea. 


Turmalinea Orihara & N. Maek., gen. nov. — MycoBank 
MB803433 


Type species. Turmalinea persicina Orihara. 


Etymology. The Latinized name, Turmalinea, refers to tourmaline, the 
reddish or bluish coloured mineral, referring to the striking colour variety of 
the fruitbodies. 


Fruitbodies mostly less than 20 mm diam, solitary to sparse, 
subglobose to depressed-globose to reniform, rubbery, sessile 
or rarely with a short stipe at the base, surface smooth to slightly 
felty, pale pink to pink, or white to brownish white, often turning 
blue to indigo-blue when rubbed or bruised. Rhizomorphs whit- 


ish or yellow to orange. Gleba initially white, then maturing to 
blackish brown, firm, loculate, of minute irregular chambers and 
thin tramal plates, chambers arranged somewhat radially near 
centre. Trama mostly pulvinate, rarely becoming a short stipe, 
subgelatinous. Columella absent. Peridium thin, up to 450 um 
in dried specimens, single- or two-layered, of interwoven to 
parallel, septate, filamentous hyphae. Subhymenium undevel- 
oped. Basidia clavate to cylindrical, 2—4-spored. Basidiospores 
statismosporic, ovoid to fusoid, inamyloid, non-dextrinoid, brick- 
red to dark brown at maturity, with 5—10 longitudinal, often 
branched, costal to irregularly broken ridges up to 3.5 um high, 
with or without a hilar appendage. Clamp connections absent 
in all tissues. 


Notes — Turmalinea is morphologically characterized by 

ovoid to fusoid, reddish brown to dark brown basidiospores 
with 5—10 irregularly furcate, longitudinal ridges, and a cushion- 
like, subgelatinous sterile base. Two of the four species of 
Turmalinea have pinkish fruitbodies with yellow to orange rhizo- 
morphs, which are unusual in hypogeous sequestrate fungi 
(Fig. 7a—d). The sister genus, Rossbeevera is morphologically 
different from Turmalinea in having a sterile base that often 
forms a short, reduced stipe and paler ellipsoid or fusoid to 
fusiform basidiospores that have 3—5 longitudinal ridges. Tur- 
malinea also tends to form firm, depressed-globose to reniform 
basidiomata, which are less common in Rossbeevera (Fig. 7). 
In addition, there is a considerable genetic divergence within 
the ITS2 minisatellite-like insertion between the two genera, 
although the insertions are found at the same position in the 
ITS2 region (Fig. 6). 
Ecologically, species of Turmalinea occur in broad-leaved fo- 
rests and are presumably symbiotic ectomycorrhizal associates 
of trees in the Fagaceae. Although we have not confirmed their 
ectomycorrhizal status, it is likely that they are mycorrhizal 
based on the fact that most species in Boletaceae form ecto- 
mycorrhizas (Nuhn et al. 2013). While species of Rossbeevera 
exhibit tropical to temperate distribution and often fruit in sum- 
mer, Turmalinea species have not yet been found in the tropics 
and they generally fruit during cooler weather. This might reflect 
physiological difference in hyphal growth between members 
of the two genera and, thus, might result in their distributional 
difference. 


Another similar leccinoid sequestrate genus, Chamonixia, also 
has basidiospores with 6—12 longitudinal ridges, but the basidio- 
spores are generally ellipsoid to broadly ellipsoid (Q = 1.3—1.8). 
Macroscopically, Chamonixia spp. form light, fragile fruitbodies 
whereas those in the Rossbeevera-Turmalinea lineage have a 
more rubbery texture. In addition, the sterile base of Chamo- 
nixia fruitbodies is never gelatinized, and often forms a per- 
current or branched columella. The genus Rhodactina is a 
sequestrate member of Boletaceae reported from India and 
Thailand (Pegler & Young 1989, Yang et al. 2006). The genus 
is similar to Turmalinea in that the fruitbodies stain reddish to 
purplish and the basidiospores have 5—10 longitudinal ridges 
(Fig. 9k, I), but is distinctive from Turmalinea in the violet 
brown to purplish carmine gleba and well-developed (3—5 um 
high), unbranched, acute ridges of the basidiospores (Yang 
et al. 2006; T. Orihara pers. obs.). The ATP6 sequence from 
Rhodactina incarnata is c. 93—94 % similar to ATP6 sequences 
of Turmalinea spp., suggesting that they are not members of the 
same lineage. We obtained an ITS1 sequence of the isotype 
of R. incarnata (GenBank KC552020), and the BLAST search 
suggests that it is most closely related to species of Tylopilus. 
The southern hemisphere genus, Austrogautieria, also has 
similar basidiospores with eight or more irregularly branched, 
longitudinal ridges. However, molecular phylogenetic analysis 
shows that Austrogautieria is a member of Gallaceaceae in the 
Hysterangiales (Hosaka et al. 2006), and therefore, is phylo- 
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Fig. 7 Fruitbodies of Turmalinea and Rossbeevera species. a. Turmalinea persicina (holotype); b. Turmalinea yuwanensis (holotype); c. aberrant fruitbody 
of T. persicina with a short stipe and incomplete lamellae (arrow) (KPM-NC 18008); d. Turmalinea chrysocarpa (holotype); e. Turmalinea mesomorpha subsp. 
mesomorpha (holotype); f. Turmalinea mesomorpha subsp. sordida (holotype); g. Rossbeevera paracyanea (holotype); h. Rossbeevera cryptocyanea (holo- 
type). — Scale bars = 1 
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genetically distant from Turmalinea. Morphologically, species 
of Austrogautieria are readily distinguished from Turmalinea by 
the glutinous spore mass and the presence of a well-developed 
columella in the gleba. 


Turmalinea persicina Orihara, sp. nov. — MycoBank 
MB803437; Fig. 7a, c, 8a—e 


Holotype. Japan, Kyoto Pref., Sakyo-ku, lwakura-agura-cho, Mt Amabuki, 
c. 200 m northwest from Jisso-in Temple, under Castanopsis sieboldii, 8 Dec. 
2008, T. Orihara, Orihara952 (KPM-NC 18001; isotype TNS-F-55010). 


Etymology. Latin, persicina (= peach-coloured), refers to the characteristic 
colour of the fruitbodies. 


Fruitbodies up to 20 mm, subglobose to depressed-globose to 
reniform, rubbery, sessile or rarely with a reduced stipe at the 
base, surface covered with thin, smooth, pink to pinkish white 
peridium occasionally turning weakly blue to bluish green when 
bruised, becoming blackish pink at maturity due to blackish 
colour of the inner gleba. Gleba off-white, occasionally turning 
bluish green in youth when cut, then blackish brown at maturity, 
rubbery, composed of minute, irregular locules. Sterile base 
present, mostly pulvinate or rarely becoming a short, reduced 
stipe, translucent, subgelatinous. Rhizomorphs orange-yellow, 
common. Odour somewhat sweet but unpleasant. Basidio- 
spores 13.1—20.2(—20.6) x (7.7—)7.8-10.4(—10.7) um, mean 
16.7 x 9.1 um (SD: 1.74 (length), 0.64 (width)), Q = 1.37-2.37, 
Q,,= 1.84, symmetric, ovoid to fusoid, colourless at first then 
becoming dark brown at maturity, with 6—10 irregularly longitu- 
dinal, partially branched ridges up to 2.9 um high in water, often 
containing one large oil drop inside, often with a distinct hilar ap- 
pendage up to 4.5 um long, HA/S = 0.09-0.25, HA/S „ = 0.17 
(n = 35), spore walls 0.6—1.3 um thick. Basidia evanescent, 
cylindrical, 17.5—30.3 x 6-8.2 um, mean 24 x 7.3 um (n = 10), 
colourless to fulvous, 2-spored, walls thin (< 0.8 um thick), inner 
matrix somewhat granulate. Basidioles persistent, 10—40 um 
diam, clavulate to subspherical, colourless, thin-walled (< 0.7 
um thick). Trama colourless, of compactly interwoven, sinuate, 
thin-walled (< 0.6 um thick) filamentous hyphae 2.3-9 um 
broad. Sterile base of interwoven, partially branched, sinuate, 
colourless to pale yellowish brown, septate, thin-walled (< 0.8 
um thick) filamentous hyphae 3-13.5 um broad. Peridium 
40-325 um thick, stramineous to yellow-brown under light mi- 
croscopy, composed of non-inflated, partially branched, septate 
filamentous hyphae 2.5—9.5 um broad parallel to subparallel to 
surface forming a cutis, walls 0.5—1.5 um thick. Rhizomorphs 
50—100 um broad, of a bundle of yellow-brown, straight, thin- 
walled (< 0.8 um thick) filamentous hyphae, mostly 2—6 um 
broad but sometimes inflated to 12 um broad. 

Habitat, Distribution & Season — Hypogeous under ever- 
green plants of the Fagaceae: Castanopsis sieboldii, C. cus- 
pidata, and Quercus glauca; Japan (western Honshu, Shikoku, 
Kyushu); late autumn to spring (November to April), occasionally 
early summer (June). 


Specimens examined. Japan, Ehime Pref., Kihoku-cho, Narukawa Gorge, 
in the Castanopsis and Quercus forest, 25 Oct. 2010, F. Nagao, Nagao 
10-10-25-02, KPM-NC 17744; Kyoto Pref., Kyoto-shi, Higashiyama-ku, Mt 
Kiyomizu, near Shogun-zuka, under Castanopsis cuspidata, 28 Nov. 2004, A. 
Kajiyama, Orihara225, KPM-NC 17851; same locality, 20 Nov. 2005, T. Ori- 
hara, Orihara354, KPM-NC 17857; same locality, 26 Nov. 2006, T. Orihara, 
Orihara555, KPM-NC 17853; same locality, 27 Apr. 2008, T. Orihara, Ori- 
hara775, KPM-NC 17854; same locality, parasitized by Sepedonium chryso- 
spermum, 27 Apr. 2008, T. Orihara, Orihara775b, KPM-NC 17857; same 
locality, 22 Feb. 2009, T. Orihara, Orihara955, KPM-NC 17858; Kyoto-shi, 
Sakyo-ku, Iwakura-agura-cho, Mt Amabuki, under Castanopsis sieboldii, 
7 Nov. 2005, Y. Kotera & T. Orihara, Orihara339, KPM-NC 17852; same 
locality, 8 Dec. 2008, T. Orihara, Orihara952 (holotype KPM-NC 18001; 
isotype TNS-F-55010); same locality, 10 Dec. 2011, T. Orihara & Y. Kotera, 
KPM-NC 18009; Mie Pref., Seki-machi, near Shoho-ji Temple, under C. sie- 
boldii, 2 June 2008, M. Ohkubo & T. Orihara, Orihara776, KPM-NC 17856; 


Hyogo Pref., Shiogahara, under C. cuspidata, 20 Nov. 2011, M. Ohmae, 
KPM-NC 18053; Hiroshima Pref., Hiroshima-shi, Higashi-kum Hukuda-cho, 
under Q. glauca, 23 Oct. 2010, A. Hadano, Orihara1291, KPM-NC 18002; 
same locality, parasitized by S. chrysospermum, 23 Oct. 2010, A. Hadano, 
Orihara1292, KPM-NC 18003; Nara Pref., Nara-shi, Mt Kasuga Primeval 
Forest, under C. sieboldii, 3 Nov. 2010, T. Orihara, Orihara1345, KPM-NC 
18004; same locality, 11 Dec. 2012, H. Inui, M. Inui & T. Orihara, KPM-NC 
18056; same locality, 12 Dec. 2012, H. Inui, M. Inui & T. Orihara, KPM-NC 
18057 (duplicate TNS-F-55012); same locality, 12 Dec. 2012, H. Inui, M. Inui 
& T. Orihara, KPM-NC 18058; Oita Pref., Saiki-shi, Shiroyama Park, under 
C. sieboldii, 15 Jan. 2011, Y. Sunada & T. Orihara, Orihara1363, KPM-NC 
18006; Saiki-shi, Nakayama, near Shirohachiman Shrine, under C. sieboldii, 
15 Jan. 2011, Y. Sunada & T. Orihara, Orihara1364, KPM-NC 18007 (dupli- 
cate: TNS-F-55011); Saiki-shi, Ume-oaza, Shigeoka, 1.5 km southwest from 
Sotaro Station, under C. sieboldii, 16 Jan. 2011, Y. Sunada & T. Orihara, Ori- 
hara1368, KPM-NC 18008; Yamaguchi Pref., Yamaguchi-shi, Tokuchi, near 
the Lake Ohara, under C. cuspidata and Q. serrata, 28 Nov. 2014, Y. Kawa- 
guchi, KPM-NC 23395. 


Notes — Turmalinea persicina, which is widely distributed 
throughout western Japan, is easily distinguished from other 
species by its pinkish white to pale pink fruitbodies with yellow to 
orange rhizomorphs at the base. The large, persistent, inflated 
basidioles are also characteristic of the species (Fig. 8b). Like 
other species of Turmalinea, its sterile base is mostly pulvinate, 
but we found one fruitbody with a distinct stipe and a trace of 
hymenophoral pores (Fig. 7c; KPM-NC 18008). However, DNA 
sequences from this specimen were nearly identical to other, 
more typical specimens collected nearby (i.e., KPM-NC 18007 
and KPM-NC 18008). This is an example of the high morpho- 
logical plasticity of sequestrate fungi, sometimes caused by 
only a few nucleotide substitutions in key functional genes, 
which might be expected during the evolution of sequestrate 
taxa from epigeous relatives (Thiers 1984, Bruns et al. 1989, 
Kendrick 1992, Bougher et al. 1993, Hibbett et al. 1994). 


Turmalinea yuwanensis Orihara, sp. nov. — MycoBank 
MB803438; Fig. 7b, 8f—j 


= Hymenogaster sp., in Yoshimi & Doi, Mem. Natl. Sci. Mus. (Tokyo) 22: 
36. 1989. 


Holotype. JAPAN, Kagoshima Pref., Amami-oshima Island, Yamato-son, 
southeast foot of Mt Yuwan, 29 Nov. 2008, T. Orihara, A. Okuda, M. Ohmae, 
Orihara942, KPM-NC 18011. 


Etymology. The epithet, yuwanensis, refers to the type locality of the spe- 
cies, Mt Yuwan-dake, which is located in Amami-oshima Island, Kagoshima 
Pref., Japan. 


Fruitbodies up to 20 mm, depressed-globose to reniform, rub- 
bery, sessile or a short, rudimentary stipe at the base, peridium 
surface smooth, vivid pink but white to pinkish white near the 
base, colour change not observed when rubbed or bruised, 
becoming dull pink at maturity due to a blackish colour of the 
inner gleba. Gleba white in youth, then blackish brown at ma- 
turity, rubbery, composed of minute, irregular locules. Sterile 
base present, mostly pulvinate, translucent, subgelatinous, 
colourless when immature then becoming dull yellowish brown. 
Rhizomorphs vivid yellow, often developed. Odour somewhat 
sweet but unpleasant, similar to that of Turmalinea persicina. 
Basidiospores 11.2—(11.4—)14.6(-14.5) x 7-(7.1—)8.8(—9) um, 
mean 12.9 x 7.9 um (SD: 0.86 (length), 0.44 (width)), Q = 1.41— 
1.84, Q = 1.64, ovoid to citriform to fusoid, colourless at first 
then becoming dark brown at maturity, with 5—10 irregularly 
longitudinal ridges 1.5—3.5 um high and 2-4.1 um wide in 
water, but the ridges mostly broken, becoming somewhat flat, 
large warts, a hilar appendage often present but relatively 
undeveloped, usually less than 3 um long, HA/S = 0.09-0.24, 
HA/S, = 0.16 (n = 20), occasionally with a hyaline pedicel up to 
5.5 um long, spore walls 0.7—1.2 um thick. Basidia evanescent, 
cylindro-clavate, 16.7—32.2 x 6.8—9 um, mean 24.8 x 8.1 um 
(n = 10), colourless, 2—4-spored, walls thin (< 0.8 um thick), 
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Fig.8 a-e: Turmalinea persicina. a. Peridial hyphae (KPM-NC 17744); b. hymenium with inflated basidioles (KPM-NC 17744); c. 2-spored basidium (KPM-NC 
18001, holotype); d. basidiospores (KPM-NC 18007); e. SEM image of basidiospores (holotype). — f—j: Turmalinea yuwanensis (holotype): Peridial hyphae; 
g—h. 3- and 4-spored basidia; i. basidiospores; j. SEM image of basidiospores. — Scale bars: a, f = 20 um; b—e, g—j = 10 um. 


inner matrix inconspicuous, contains several oil drops, ster- 
igmata 4.5-7 um long. Basidioles persistent, 8—20 um diam, 
clavulate to subspherical to vesiculate, colourless, thin-walled 
(< 0.8 um thick). Trama white colourless, of compactly inter- 
woven, thin-walled (< 0.8 um thick) filamentous hyphae 2—9 
um broad. Sterile base of compactly interwoven, thin-walled 
(< 0.8 um broad), colourless to pale yellowish brown, septate 
filamentous hyphae 3—10.5 um broad. Peridium 150—400 um 
thick, in the outer part stramineous to yellow-brown under light 
microscopy, in the inner part almost colourless, composed of 
interwoven to subparallel, partially branched, non-inflated, 
septate filamentous hyphae 3—10(—13.5) um broad forming a 
cutis, walls up to 1um thick. 

Habitat, Distribution & Season — Hypogeous to subepige- 
ous under or on rotten wood of Castanopsis sieboldii subsp. 
lutchuensis; So far known from Amami-oshima, Okinawa and 
Tokashiki Islands, the Ryukyu Archipelago, Japan; late autumn 
(November). 


Specimens examined. Japan, Kagoshima Pref., Amami-oshima Island, 
Naze-shi, Kinsakubaru National Forests, 19 Nov. 1988, Y. Doi, ‘Hymenogaster 
sp.’, Yoshimi7393, TNS-F-183213; Amami-oshima Island, Yamato-son, 
southeast foot of Mt Yuwan, under Castanopsis sieboldii subsp. /utchuensis, 
17 Nov. 2007, A. Hadano, T. Orihara, M. Ohkubo, Orihara758, KPM-NC 
18010 (duplicate: TNS-F-55013); same locality, 29 Nov. 2008, T. Orihara, 
A. Okuda, M. Ohmae, Orihara942, KPM-NC 18011 (holotype); Okinawa Pref., 
Tokashiki-jima Island, Tokashik-son, Tokashiki, under Castanopsis sieboldii 


subsp. /utchuensis, 14 Apr. 2013, T. Orihara, KPM-NC 23377; Okinawa 
Island, Kunigami-son, along the path to Mt Yonaha, under Castanopsis 
sieboldii subsp. lutchuensis, 10 Nov. 2013, M. Moriguchi, KPM-NC 23390; 
same locality, 5 Feb. 2015, T. Orihara, KPM-NC 23998; same locality, 5 Feb. 
2015, T. Orihara, KPM-NC 24000. 


Notes — The vivid pink peridial surface and yellow rhizo- 
morphs of Turmalinea yuwanensis are striking features and are 
easily recognized in the field. Microscopically, T. yuwanensis is 
unique in that its basidiospores have distinct but discontinuous 
longitudinal ridges that partially become large, irregular warts. 
Asimilar and closely related species, T. persicina, is distinguished 
from T. yuwanensis by the paler peridium, larger basidiospores 
with relatively continuous ridges, and its slender, 2-spored 
basidia. The Chinese species, T. chrysocarpa is also phyloge- 
netically close to T. yuwanensis, but is readily distinguished by 
the yellow peridium and larger basidiospores (Fig. 8d). 


This species has been collected from only a few sites in the 
central Ryukyu Islands, which lie c. 300—600 km to the south- 
west of mainland Japan. A truffle-like fruitbody with orange 
rhizomorphs reported by Yoshimi & Doi (1989) from Amami- 
oshima island as ‘Hymenogaster sp.’ was confirmed as an im- 
mature T. yuwanensis. Given the unique phylogenetic position 
and rarity of T. yuwanensis, the conservation of its habitat is 
of primary importance in terms of phylogenetic diversity (Faith 
1992, Faith & Baker 2006). 
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Turmalinea chrysocarpa Orihara & Z.W. Ge, sp. nov. — Myco- 
Bank MB803441; Fig. 7d, 9a—c 


Holotype. Cina, Yunnan Prov., Manhao, Gejiu, alt. 860 m, in a broad- 
leaved forest, 24 Sept. 2011, Z.W. Ge 3098 (HKAS70601; isotype KPM-NC 
18068). 


Etymology. The Latin, chryso- (golden) and -carpa (fruitbody), refer to the 
colour of the fruitbody surface of the species, which is unusual in Turmalinea. 


Fruitbody solitary, 1.5 cm diam, depressed-globose, surface 
smooth to slightly floccose, yellow to light orange, stipe reduced, 
rudimentary, with orange rhizomorphs at the base. Gleba choco- 
late brown to blackish brown at maturity, composed of minute, 
irregular locules. Sterile base present, more or less pulvinate, 
translucent, subgelatinous. Odour strong and sweet. Basidio- 
spores 13.5—(13.8—)21.5(—23.5) x (6.9—)7.3—9.7(-10.8) um, 
mean 17.5 x 8.5 um (SD: 1.97 (length), 0.61 (width)), Q = 1.66— 
2.71, Q, = 2.06, symmetric, citriform to fusiform, dark reddish 
brown at maturity, with 6—9 irregularly longitudinal, often furcate 
ridges 1—2.2 um high in water, with a distinct hilar appendage 
up to 6 um long, HA/S = 0.11-0.31, HA/S „ = 0.19, spore walls 
0.7—1.2 um thick. Basidia evanescent, cylindrical but mostly 
collapsed, 20.9 x 3.2 um, 2-spored, colourless, sterigmata up 
to 7.5 um long. Trama colourless, of moderately or somewhat 
loosely interwoven, partially sinuate, septate filamentous hy- 
phae 2—10 um broad, walls up to 1.2 um thick. Peridium 150— 
330 um thick, composed of two layers: outermost layer (i.e., 
peridiopellis) 60—150 um thick, of parallel to subparallel, 
straight, septate, colourless, thin-walled (< 0.8 um), non-inflated 
filamentous hyphae 2.8—10(—13) um broad, only the surface 
stained fulvous; inner layer (i.e., subpellis) 75—250 um thick, 
colourless, of compactly interwoven, sinuate, septate, narrower 
(2.5—8 um), thin-walled filamentous hyphae. 


Notes — Turmalinea chrysocarpa, which is reported from 
Southwest China, is unique in having a yellow to light orange 
peridium and an inner peridial layer composed of compact, 
interwoven hyphae. Phylogenetic analyses support that T. chry- 
socarpa is a distinct lineage from the other Turmalinea spp. 
from Japan (Fig. 3-6). This species is readily distinguished 
from other members of the genus by the colour of the peridium 
and the dimension of basidiospores. Unfortunately, only one 
specimen of T. chrysocarpa has so far been collected. Although 
this may be due to the rarity of the species, it may also be 
due to lack of data on the diversity, phenology and ecological 
requirements of sequestrate fungi in China. Despite our lack 
of knowledge about this species, its characteristic morphology 
and phylogenetic distinctness are sufficient to describe it as a 
new Turmalinea species. 


Turmalinea mesomorpha Orihara, sp. nov. — MycoBank 
MB803439; Fig. 7e, 9d—g 


Holotype. Japan. Aomori Pref., Towada-shi, Ohoronai, under Fagus cre- 
nata, 12 Oct. 2009, T. Sasa & T. Orihara, Orihara1080, KPM-NC 18015. 


Etymology. The Latin, meso- (middle or intermediate) and -morpha (mor- 
phology), signifying that the species has intermediate morphology between 
Rossbeevera spp. and other members of Turmalinea. 


Subspecies autonym. Turmalinea mesomorpha subsp. mesomorpha (the 
automatically generated name by the proposal of T. mesomorpha subsp. 
sordida, which is to be described below). 


Fruitbodies up to 24 mm, subglobose to depressed-globose 
to reniform, rubbery, sessile or with a short, reduced stipe less 
than 5 mm long at the base, surface somewhat felty or partially 
wrinkled, white to off-white, immediately turning indigo to pur- 
plish blue and finally blackish brown when touched or bruised. 
Gleba off-white to whitish yellow when immature, then chocolate 
brown, and finally becoming blackish brown at maturity, rubbery, 
composed of minute, irregular locules. Sterile base present, 


mostly pulvinate but rarely columella-like, translucent, subge- 
latinous, colourless to blue-grey when immature then becoming 
dull yellowish brown at maturity, surrounded by whitish hyphal 
veins irregularly radiating from the central sterile base and 
connected to inner portion of the peridial context at the bottom. 
Rhizomorphs white to off-white, common, showing the same 
pattern of discolouration as the fruitbody surface. Odour some- 
what sweet but unpleasant. Basidiospores 12.5—16.4(-16.9) x 
(7.8—)7.9-10.4 um, mean 14.4 x 9.1 um (SD: 0.97 (length), 0.61 
(width)), Q = 1.38-1.92, Q, = 1.58, symmetric, ovoid to citriform 
to fusoid, colourless at first then becoming dark reddish brown 
at maturity, with 5—10 irregularly longitudinal, costal or often 
branched ridges up to 2.6 um high in water, a hilar appendage 
reduced or indistinct, occasionally with a hyaline pedicel up to 
6 um long, spore walls 0.7—1.6 um thick. Basidia evanescent, 
clavate to cylindro-clavate, 25.5—43.8 x 6.4—12.2 um, mean 
33.1 x 8.8 um (n = 15), 2—4-spored, colourless or with a yellow- 
ish brown tint, sterigmata 6—10.5 um long, walls thin (< 0.8 um 
thick), inner matrix somewhat granulate. Basidioles cylindrical 
to cylindro-clavate, 10—18 x 6—9 um, colourless or pigmented 
fulvous, walls thin (< 0.7 um thick). Trama colourless or pig- 
mented pale yellowish brown in places, of subparallel to loosely 
interwoven, straight, partially branched, thin-walled (< 0.8 um 
thick) filamentous hyphae 2—10 um broad. Sterile base of 
loosely interwoven, mostly sinuous, colourless, septate, thin- 
walled (> 0.8 um thick) filamentous hyphae 3.5—11 um broad, 
the surrounding whitish layer composed of loosely interwoven, 
more or less frequently branched, thin-walled (< 0.8 um thick), 
filamentous hyphae 2—8 um broad that are less sinuate than 
those comprising the sterile base. Peridium 50—150 um thick, 
yellow-brown under light microscopy, composed of repent, 
partially branched and inflated, septate filamentous hyphae 
2-11 um broad interwoven or subparallel to surface, walls up 
to 1 um thick. 

Habitat, Distribution & Season — Hypogeous to subepige- 
ous under Fagus crenata; Japan (northern Honshu); autumn 
(September to October). 


Specimens examined. Japan, Aomori Pref., Towada-shi, Ohoronai, under 
Fagus crenata (the holotype locality), 21 Sept. 2008, T. Sasa, Orihara874, 
KPM-NC 18012 (duplicate: TNS-F-55014); same locality, 12 Oct. 2009, T. Sasa 
& T. Orihara, Orihara1080, KPM-NC 18015 (holotype); Towada-shi, Tsuta- 
onsen, under F. crenata, 27 July 2012, T. Muroi, KPM-NC 18059; Iwate Pref., 
Hachimantai-shi, Appi-kohgen, under F. crenata, 11 Oct. 2009, M. Taniguchi, 
Orihara1074, KPM-NC 18013; same locality, 12 Oct. 2009, T. Sasa & T. Ori- 
hara, Orihara1075, KPM-NC 18014. 


Notes — Turmalinea mesomorpha, the earliest diverging 
species within Turmalinea, shares characteristics of both Tur- 
malinea and Rossbeevera. Like most Rossbeevera species, 
the fruitbody has a whitish peridium that quickly turns indigo 
to purplish blue. On the other hand, most other generic level 
characters such as the cushion-like sterile base and the basi- 
diospore morphology are typical of Turmalinea. The morpho- 
logy of T. mesomorpha in combination with its position in the 
phylogeny suggests that the pinkish to yellow colours of the 
other Turmalinea species probably evolved after the divergence 
of T. mesomorpha. Compared to other Turmalinea species, 
the sterile base tends to become amorphous and surrounded 
with whitish hyphal veins. Interestingly, one fruitbody (KPM-NC 
18059) had a columnar sterile base, which is more common in 
Rossbeevera. These peridial and sterile base characters mor- 
phologically support that the species is the earliest diverging 
lineage within the genus Turmalinea. Moreover, the aberrant 
sterile base morphologies found in T. mesomorpha and T. per- 
sicina suggest that these reversions have occurred somewhat 
frequently within the genus. 


Turmalinea mesomorpha is microscopically different from T. per- 
sicina and T. yuwanensis in its larger basidia and medium-sized 
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Fig.9 a-c: Turmalinea chrysocarpa (HKAS 70601, holotype). a. Peridium (arrows: pp = outer layer (peridiopellis); il = inner layer); b. SEM image of basidio- 
spores; c. basidiospores. — d—g: Turmalinea mesomorpha subsp. mesomorpha. d. Peridial hyphae (KPM-NC 18014); e. basidia (KPM-NC 18014); f. SEM 
image of basidiospores (KPM-NC 18015, holotype); g. basidiospores (holotype). — h-j: Turmalinea mesomorpha subsp. sordida. h. Peridial hyphae (KPM- 
NC 17743); i. SEM image of basidiospores (KPM-NC 18016, holotype); j. basidiospores (holotype). — k—I: Basidiospores of Rhodactina spp. k. Rhodactina 
himalayensis; |. Rhodactina incarnata. — Scale bars: a = 50 um; b, c, e—g, i—l = 10 um; d, h = 20 um. 


basidiospores with relatively continuous longitudinal ridges. In 
contrast to T. persicina and T. yuwanensis, which occur in the 
evergreen Quercus/Castanopsis forest, T. mesomorpha has 
only been collected under the deciduous montane tree Fagus 
crenata. This suggests not only geographical but also ecologi- 
cal isolation from the other two Japanese Turmalinea species. 


Turmalinea mesomorpha subsp. sordida Orihara, subsp. 
nov. — MycoBank MB803440; Fig. 7f, 9h-j 


Holotype. JAPAN, Ehime Pref., Saijo-shi, Mt Ishizuchi, under F. crenata, 
10 Oct. 2010, H. Ikeda, Orihara1280, KPM-NC 18016. 


Etymology. The Latin, sordida (= dirty, foul), refers to the characteristic 
colour of the peridial surface of the subspecies. 


Diagnosis. This subspecies is distinguished from Turmalinea mesomorpha 
subsp. mesomorpha by the greyish brown surface of fruitbodies that 
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turns deep blue to purple and finally becoming blackish brown when 
touched or bruised, and thicker peridium up to 450 um thick composed of 
filamentous hyphae 3—15 um broad. Basidiospores shape and dimensions 
not significantly different from that of T. mesomorpha subsp. mesomorpha, 
11.7—(12—)15.6 (-17) x (6.8—)7.3—9.6 um, mean 13.7 x 8.5 um (SD: 0.96 
(length), 0.58 (width)), Q = 1.41-1.94, Q, = 1.62. 


Habitat, Distribution & Season — Hypogeous to subepige- 
ous under Fagus crenata; Japan (Shikoku); autumn (September 
to October). 


Specimens examined. Japan, Ehime Pref., Saijo-shi, Mt Ishizuchi, under 
F. crenata, 10 Oct. 2010, H. Ikeda, Orihara1280, KPM-NC 18016 (holotype); 
Matsuyama-shi, Mt Takanawa, under F. crenata, 4 Nov. 2007, F. Nagao, 
KPM-NC 17743. 


Notes — The T. mesomorpha specimens collected from 
Shikoku, one of the four main islands of Japan, showed slight 
macroscopic differences from the other collections from the 
northern part of Japan, including the holotype specimen of 
subsp. mesomorpha. However, we did not find any significant 
microscopic differences between the two groups, and they 
share a similar ecological habitat (i.e., they are all found in as- 
sociation with Japanese beeches). The phylogenetic analyses 
conducted in this study showed slight divergence between the 
two lineages, suggesting their relatively recent diversification 
compared to the other species-level lineages in the genus. We 
therefore recognize the Shikoku lineage as subspecies sordida 
based on morphological and molecular data that are currently 
available. 


Aoki (1978) informally reported a specimen from Tokyo that 
more or less fit the characters of T. mesomorpha subsp. sor- 
dida as ‘Gautieria sp.’ However, we were unable to locate this 
specimen, so we could not directly compare it to our recently 
collected material. Since Tokyo is geographically intermediate 
between the present habitats of each subspecies, the record 
or additional collections in that area should provide important 
knowledge about the biogeography of the species in the future. 


Rossbeevera T. Lebel & Orihara, in Lebel et al., Fung. Diversity 
52: 54 + 73. 2012, ‘Rosbeeva’. 


Type species. Rossbeevera pachydermis T. Lebel, Fung. Diversity 52: 64. 
2012, ‘Rosbeeva pachyderma’. 


Notes — The precise description of Rossbeevera is pre- 
sented in Lebel et al. (2012a) under the variant ‘Rosbeeva’. 
Since the correction of the type species name is omitted in the 
erratum by Lebel et al. (2012b), we herein specify its correct 
name. 


Rossbeevera paracyanea Orihara, sp. nov. — MycoBank 
MB803442; Fig. 7g, 10a—e 


Holotype. JAPAN, Nara Pref., Nara-shi, Nara Park, under Quercus gilva, 
26 Oct. 2008, M. Kawai, Orihara908, KPM-NC 17847. 


Etymology. The epithet (Greek, para- = para- or beside and cyanea = ul- 
tramarine blue) expressing the morphological characteristics, habitat and 
distribution of the species similar to those of Rossbeevera eucyanea: both 
of the species have strongly cyanescent fruitbodies and occur usually under 
Castanopsis and evergreen Quercus trees in western Japan. 


Fruitbodies solitary or in small clusters, up to 27 mm diam, sub- 
globose to depressed-globose, soft, peridial surface smooth, 
white to greyish at first, becoming dark blue-grey to dark grey 
and occasionally cracked at maturity to expose inner gleba and 
turning deep blue when touched or bruised. Gleba off-white 
when immature, dark reddish brown at maturity, non-gelatinous, 
minutely and irregularly loculate, turning indigo-blue in some 
portions when cut and exposed to air. Sterile base present but 
reduced, non-gelatinous, strongly cyanescent, often forming 


a short stipe. Columella absent. Odour somewhat sweet but 
unpleasant. Basidiospores (13.3—)14-—19.3(—20.4) x (6.7-) 
6.9-9.2(—9.3) um, mean 16.7 x 8 um (SD: 1.32 (length), 0.57 
(width)), Q = 1.7-2.5, Q, = 2.1, symmetric, fusoid to fusiform, 
inamyloid, nondextrinoid, colourless at first then becoming red- 
dish brown at maturity, with 4—5 or rarely 3 smooth, longitudinal 
ridges up to 2.5 um high in water, walls 0.5—1.2 um thick, with 
a developed hilar appendage 2.1—4.5(—4.7) um long (mean 
3.3 um) at the base, HA/S = 0.13-0.25, HA/S „= 0.2, walls at 
the tips thinner than on the sides. Basidia 13.8—21 x 6—11 um 
(mean = 18 x 8.1 um; n = 10), clavulate to doliform, evanescent, 
colourless or pigmented with yellow-brown, mostly 2-spored 
but rarely 1-, 3- or 4-spored. Hymenium developed when im- 
mature but collapsed at maturity, colourless; basidioles clavate 
to clavulate to doliform almost the same size as or smaller than 
basidia. Subhymenium not developed, composed of single- or 
two-layered, colourless, isodiametric to clavulate cells 4—15 um 
diam. Trama of parallel to subparallel, non-inflated, colourless, 
thin-walled (< 0.6 um thick) filamentous hyphae 2.5-5.5 um 
broad. Sterile base composed of compact, interwoven, septate, 
inflated, thin-walled (c. 1 um thick) filamentous hyphae 2—18 
um broad. Peridium thin, composed of 2 layers: outer layer 
50—160 um thick, of compact, interwoven, septate, pale brown 
to brown filamentous hyphae 2—9.5 um swollen in some portion 
to 15 um broad, walls up to 1.2 um thick; inner layer 15—55 um 
thick, tan to reddish brown, of subparallel or tangled, sinuate, 
non-inflated filamentous hyphae 1.5—5 um broad, walls 0.5—1.3 
um thick. Clamp connections absent in all tissues. 

Habitat, Distribution & Season — Hypogeous to subepige- 
ous under Quercus gilva and Castanopsis cuspidata; Japan 
(western Honshu); summer to early winter (July to December). 


Specimens examined. Japan, Nara Pref., Nara-shi, Nara Park, under 
Quercus gilva, 26 Oct. 2008, M. Kawai, Orihara908, KPM-NC 17847 (holo- 
type); Nara-shi, Mt Kasuga Primeval Forest, under Quercus gilva, 24 Dec. 
2011, T. Orihara, KPM-NC 18022; same loclity and date, T. Orihara, KPM-NC 
18023; same locality, 12 Dec. 2012, T. Orihara, H. Inui, M. Inui, KPM-NC 
18087; same locality, 16 Nov. 2013, M. Ohmae, KPM-NC 23940; same 
locality, 19 July 2014, M. Ohmae, KPM-NC 23953; same locality and date, 
M. Ohmae, KPM-NC 23954; Osaka Pref., Minoh-shi, Mt Minoh, 5 July 2010, 
M. Ohmae, Orihara1178, KPM-NC 17848; same locality, 7 Aug. 2011, M. Ohmae, 
KPM-NC 18027; same locality, 21 July 2014, R. Nakano, KPM-NC 23955; 
same locality, 15 July 2014, M. Ohmae, KPM-NC 23956. 


Notes — Fruitbodies of Rossbeevera paracyanea are macro- 
scopically similar to those of another Japanese species, R. eu- 
cyanea, especially when young. Its habitat and size of basidio- 
spores also overlap with those of R. eucyanea. However, R. para- 
cyanea is morphologically distinguished from R. eucyanea by 
its two-layered peridium comprised of the outer layer of partially 
inflated, filamentous hyphae and the thin, inner layer of some- 
what pigmented, sinuate filamentous hyphae. The peridium 
tends to peel off at maturity due to the fact that the hyphae of 
the inner peridial layer are more or less parallel to the surface 
of the gleba. The surface of R. paracyanea fruitbodies was 
generally greyish or blue-grey and this helps to distinguish 
the species from R. eucyanea, which has a whitish fruitbody. 


The three-locus nuclear phylogeny showed distinct divergence 
within R. paracyanea (i.e., lineages 1 and 2), but the mitochon- 
drial phylogeny showed a conflicting relationship (Fig. 5). In the 
R. paracyanea lineage 2 specimen from Osaka Pref. (KPM-NC 
17848), some of the basidia have unusually long sterigmata (up 
to 20 um long) and the context of the outer layer of the peridium 
was almost colourless. The latter observation, however, could 
be due to the immature state of the specimens. Thus, we leave 
the species as a single taxon and treat the collections from 
Osaka (KPM-NC 17848 and KPM-NC 18027) as an infraspecific 
cryptic lineage until sufficient collections become available for 
additional phylogenetic and morphological comparisons. 
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Fig. 10 a-e: Rossbeevera paracyanea (a-d: KPM-NC 17847, holotype; e: KPM-NC 23940). a. Peridium (upper right). The inner peridial layer is indicated 
by an arrow; b. basidia and basidiospores; c. mature basidiospores; d. basidiospores mounted in lacto-glycerol after pre-soaking in 3 % KOH; e. SEM image 
of basidiospores. — f—i: Rossbeevera cryptocyanea (KPM-NC 23928, holotype). f. Peridium; g. basidiospores; h. basidiospores extending from a 4-spored 
basidium; i. SEM image of basidiospores. — Scale bars: a, f = 50 um; b—e, g—i = 10 um. 


Rossbeevera cryptocyanea Orihara — MycoBank MB811526; 
Fig. 7h, 10f-i 


Holotype. Japan, Kakeroma Isl., Setouchi-machi, Akitoku, under Castanop- 
sis sieboldii subsp. lutchuensis, 28 June 2014, T. Orihara, KPM-NC 23928. 


Etymology. The epithet (Greek, crypto- = crypto- or hidden, and cyanea = 
ultramarine blue) expressing its close phylogenetic relationship and extreme 
morphological similarity to R. eucyanea. 


Diagnosis. Morphologically quite similar to Rossbeevera eucyanea, but dif- 
fers in smaller and narrower basidiospores with longer hilar appendices. 


Fruitbodies solitary or in small clusters, up to 22 mm, subglo- 
bose to depressed-globose to reniform, soft, surface smooth 
to somewhat felty, white or more rarely greyish white, quickly 
turning deep blue when touched or exposed to air. Gleba off- 
white when immature, then becoming brown, non-gelatinous, 
minutely and irregularly loculate, quickly turning indigo-blue 
when cut and exposed to air. Sterile base present but reduced, 
non-gelatinous, cyanescent, forming a short stipe. Columella 
absent. Odour somewhat sweet but unpleasant. Basidiospores 
13.4-—(13.5—-)18.3(-18.4) x 5.8(-7.2)—7.3 um, mean 15.8 x 
6.5 um (SD: 1.2 (length), 0.36 (width)), Q = 1.9-3.0, Q, = 2.4, 


fusoid to fusiform, colourless at first then becoming dark brown, 
with 3—5(—6) longitudinal ridges up to 2 um high in water, walls 
0.5-1.2 um thick, with large, developed hilar appendix 2.4—5.9 
um long (mean 3.7 um) at the base, HA/S = 0.17—0.37, HA/S „ 
= 0.23. Spore walls at the tips tend to be thinner than on the 
sides. Basidia 12.5—21 x 7-11.5 um, clavulate to clavate, 
colourless, 2- or 4-spored. Trama composed of interwoven, 
septate, colourless, thin-walled, filamentous hyphae 3—8 um 
broad. Peridium 50-180 um thick in dried specimens, of some- 
what loosely woven, septate, partially branched or swollen, 
thin-walled (< 0.5 um) filamentous hyphae 3—15 pm, subparal- 
lel to surface, surface often stained blue. Clamp connections 
absent in all tissues. 

Habitat, Distribution & Season — Hypogeous to subepigeous 
under Castanopsis sieboldii subsp. /utchuensis and Quercus 
glauca; Japan (Kyushu and the Ryukyu Archipelago); early 
summer and late autumn (June to November). 


Specimens examined. Japan, Oita Pref., Usuki-shi, Shimonoe, under Quer- 
cus glauca, 20 Nov. 2005, Y. Sunada, Orihara1351, KPM-NC 17846; Kago- 
shima Pref., Amami-oshima Isl., Yamato-son, north-eastern foot of Mt Yu- 
wandake, under Castanopsis sieboldii subsp. lutchuensis, 17 Nov. 2007, 
A. Hadano & T. Orihara, Orihara760, KPM-NC 17843; Kakeroma Isl., Setouchi- 
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machi, Akitoku, under Castanopsis sieboldii subsp. Jutchuensis, 28 June 
2014; T. Orihara, KPM-NC 23928 (holotype); Okinawa Pref., Ishigaki Isl., 
Ishigaki-shi, near Omoto tunnel, under Castanopsis sieboldii subsp. /utchu- 
ensis, 11 Oct. 2013, M. Suyama, KPM-NC 23387. 


Notes — Macroscopically, R. cryptocyanea is indistinguish- 
able from R. eucyanea: they both have white fruitbodies immedi- 
ately discolouring indigo-blue to ultramarine when touched and 
share the same habitat in the southern part of Japan (i.e., they 
both occur under evergreen trees of the Fagaceae). The only 
reliable diagnostic character is the basidiospore morphology. 
Basidiospores of R. cryptocyanea are narrower than that of 
R. eucyanea (Q = 1.9-3 vs 1.6-2.2; Q. = 2.4 vs 1.9), and the 
hilar appendix of R. cryptocyanea basidiospores is much longer 
than that of R. eucyanea. This results in a distinct difference 
between their HA/S values and, hence, their basidiospore 
proportion: HA/S of R. eucyanea is significantly smaller than 
that of R. cryptocyanea (i.e., HA/S = 0.11-0.18, HA/S „= 0.14 
in R. eucyanea; n = 25). 


DISCUSSION 


Phylogeny, evolution, phylogeography and morphology of 
Rossbeevera and Turmalinea 


In this study, we took several different phylogenetic approaches, 
including nuclear and mitochondrial multi-locus analyses, to 
clarify evolutionary relationships of Rossbeevera and the allied 
taxa in the leccinoid clade. The nuclear three-locus phylogeny 
based on the combined dataset confirmed that sequestrate 
fungi fall into several distinct lineages within the leccinoid 
clade and this phylogeny also elucidated the placement of the 
new generic sequestrate lineage, Turmalinea (Fig. 3). Loss of 
forcible spore discharge (i.e., ballistospory) on the process of 
evolution of sequestrate morphology is thought to be irreversi- 
ble (Hibbett 2004, Wilson et al. 2011). The major morphologi- 
cal characteristics of the sequestrate genera Chamonixia and 
Octaviania, such as basidiospores and glebal structure, are 
quite different from each other and from Rossbeevera and Tur- 
malinea (Lebel et al. 2012a, Orihara et al. 2012b). The nuclear 
multi-locus combined phylogeny shows that the sequestrate 
genera in the leccinoid clade are not monophyletic (Fig. 3). 
In contrast, the overall morphology between the two boletoid 
genera, Leccinum and Leccinellum, is similar (Den Bakker & 
Noordeloos 2005). Unfortunately, the phylogenetic placement 
of Leccinellum in the leccinoid clade has not been settled and 
the placement of Chamonixia as the earliest diverging lineage 
in the clade is only moderately supported (Fig. 3). Nonethe- 
less, given that the phylogenetic analysis is consistent with the 
morphological differences among the sequestrate genera and 
the fact that the dominant fruitbody morphology in Boletaceae 
is the boletoid form (Binder & Hibbett 2006, Nuhn et al. 2013), 
these results suggest that the Chamonixia, Octaviania, and 
Rosbeevera/Turmalinea lineages evolved independently from 
different boletoid ancestors in the leccinoid clade. 


Unfortunately, the multi-locus phylogenies in this study did not 
clarify the sister clade to the leccinoid clade (Fig. 2, 3) so we 
were unable to definitively determine the sequence of morpho- 
logical changes that may have occurred during the evolution 
of this group. Recently, Wu et al. (2014) proposed several 
subfamilies across the Boletaceae based on nuclear multi- 
gene (nLSU, EF-1a, RPB1, and RPB2) analyses. Among the 
proposed subfamilies was the Leccinoideae, which includes 
the genera of the leccinoid clade. Our phylogenetic analyses 
did not recover Leccinoideae but this is likely due to the dif- 
ference of loci included in the combined dataset. Instead, two 
genera in the Leccinoideae (i.e., Borofutus and Spongiforma) 
were resolved in a monophyletic clade that included Por- 
phyrellus porphyrosporus, which was treated as a member of 


the Boletoideae in Wu et al. (2014). This suggests a need to 
re-evaluate the phylogenetic entity of the Leccinoideae based 
on different loci and more inclusive boletoid and sequestrate 
taxa. The addition of other DNA loci such as RPB1 and RPB2 
will also help to resolve the generic relationships within the 
leccinoid clade more definitively. 


The phylogenetic analyses also revealed a number of new 
species-level sequestrate lineages within Rossbeevera and 
Turmalinea with robust statistical support. Accordingly, we 
propose two new Rossbeevera species and four new species 
and one new subspecies of Turmalinea. The description of 
these new taxa almost doubles the number of species within 
the Rossbeevera-Turmalinea lineage, suggesting that this 
group is likely diverse and ecologically important in East Asia, 
SoutheastAsia, and Australasia despite the fact that this lineage 
was only recently recognized (Lebel et al. 2012a, b). So far, Tur- 
malinea accommodates only a few species (4 spp.) compared 
to Rossbeevera (10 spp.), and Octaviania (c. 30 spp.), but the 
long branches in the multi-gene phylogenies (Fig. 2, 3) imply 
that there are more undescribed species yet to be discovered. 
It is also possible, however, that the long branches are the 
result of extinction of ancestral species. Collecting surveys of 
sequestrate fungi across extensive areas in East Asia will give 
further insight into the diversity of Turmalinea. 


While Turmalinea spp. are presently known only from East Asia 
(e.g., China and Japan), species of Rossbeevera have also 
been reported from Australasia and Southeast Asia (Lebel et 
al. 2012a, Orihara et al. 2012a, Vernes & Lebel 2011). Orihara 
et al. (2012a) and this study (Fig. 2, 3) agree that the Chinese 
species R. yunnanensis forms the basal branch within Ross- 
beevera. Based on these results we postulate that the genera 
Turmalinea and Rossbeevera may have originated in China or 
SE Asia and then migrated independently into Japan (including 
the Ryukyu Archipelago). 


Similarly, all the phylogenies shown in this study suggest that 
the Australasian Rossbeevera species diversified relatively 
recently and that they probably migrated from Asia only once. 
Since sequestrate fungi are thought to rely on animal myco- 
phagy for spore dispersal (Claridge & May 1994, Fogel & Trappe 
1978, Maser 1978, Maser et al. 2008, Vernes & Dunn 2009), 
long-distance spore dispersal over ocean barriers would seem 
challenging. However, a biogeographic study by Hosaka et al. 
(2008) demonstrated that intercontinental, overseas dispersal 
has occurred separately in at least three truffle-like genera 
in the Hysterangiales, suggesting that overseas dispersal of 
truffles may be more common than previously thought. Field 
observations of truffle consumption by mammals indicate that 
Rossbeevera species are highly sought after as food. For 
example, Rossbeevera spores are among the most common 
spore types found in a diverse range of small Australian mam- 
mals (Lebel unpubl. data) and Rossbeevera spores have also 
been documented in mycophagous mammals from Papua New 
Guinea (Vernes & Lebel 2011). These observations suggest 
that intercontinental migration of Rossbeevera from Asia to 
Australasia may have been facilitated by mycophagous animals. 


The sister sequestrate genera, Turmalinea and Rossbeevera, 
have several characteristics in common, including basidio- 
spores with longitudinal ridges and fruitbodies that turn bluish 
upon bruising or exposure to air. However, Turmalinea species 
have firm, rubbery fruitbodies with a pulvinate sterile base and 
basidiospores with 5—10 longitudinal ridges whereas Rossbee- 
vera species generally have softer fruitbodies with a reduced 
stipe-columella and basidiospores with 3—5 longitudinal ridges. 
In addition, three of the four Turmalinea species have brightly 
coloured, pinkish or yellow fruitbodies with yellow to orange 
rhizomorphs, whereas all Rossbeevera species have whitish to 
greyish white fruitbodies (Fig. 7). The one exceptional species 
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within Turmalinea, the earliest diverging species T. mesomor- 
pha, forms whitish fruitbodies with bluish discolouration similar 
to species of Rossbeevera. This parsimoniously suggests that 
the ancestor of the two genera had a whitish sequestrate fruit- 
body with bluish discolouration and longitudinally ridged basidio- 
spores as synapomorphies. 


Although significantly supported in the phylogenetic analyses, 
it is difficult to infer close evolutionary relationships between 
bolete genera (i.e., Leccinum and Leccinellum) and sequestrate 
genera (i.e., Rossbeevera, Turmalinea, Octaviania, and Chamo- 
nixia) in the leccinoid clade based on morphology. For exam- 
ple, basidiospores of Leccinellum and Leccinum are smooth 
and colourless to pale brown, whereas those of the sequestrate 
taxa in Rossbeevera, Turmalinea, Chamonixia, and Octavia- 
nia are distinctly ornamented and brown to blackish brown. 
In addition, the epigeous boletoid taxa in the leccinoid clade 
have cheilo- and pleurocystidia in their tissues (Den Bakker & 
Noordeloos 2005) while the sequestrate taxa generally lack 
cystidia. The only distinct character that unites these genera 
is the reddish or bluish discolouration of fruitbodies, but this 
feature can also be found widely across both sequestrate 
and epigeous boletoid taxa in the Boletaceae. Therefore, we 
suggest that the evolution of dark pigmentation and enhanced 
ornamentation of basidiospores as well as the loss of cystidia 
are likely features that are selected for during the transition 
from boletoid to sequestrate fruiting form in the leccinoid clade. 
Future analyses using phylogenetic analysis in conjunction with 
molecular dating should examine the relationships between the 
micro- and macroscopic morphologies across the Boletales and 
particularly within poorly resolved groups within Boletaceae. 


Interpreting evolutionary history in light of topological 
discord among phylogenies 


The single-locus and combined phylogenies showed consider- 
able topological discord at various degrees of divergence within 
Rossbeevera. The three single-locus nuclear trees exhibit dif- 
ferent patterns of divergence among Japanese Rossbeevera 
spp. (Fig. 1). This topological incongruence can be explained 
by deep coalescence or incomplete lineage sorting (ILS) (Mad- 
dison 1997). If this is true then sequences of ITS, nLSU, and 
EF-1a were divergent in the ancestral Rossbeevera species 
before speciation events that led to the extant Japanese and 
Australasian Rossbeevera spp. The *BEAST analysis estimates 
a Bayesian species tree that accounts for effects of ILS under 
the multispecies coalescent model (Heled & Drummond 2010). 
Most conflicting topologies in the three nuclear gene trees 
were not satisfactorily reconciled in the *BEAST species tree 
based on the three loci, but the sister relationship between 
R. eucyanea and R. cryptocyanea was strongly supported by 
the *BEAST species tree (Fig. 2). These results suggest that 
gene trees based on other loci may mislead interpretations of 
the relationship between the two species. In addition, although 
most species-level relationships in Rossbeevera were strongly 
supported in the combined phylogeny (Fig. 3), the species 
tree analysis suggests that these internal relationships must 
be viewed critically in light of the incongruence between the 
various datasets. 


We also compared the phylogeny based on three nuclear loci 
(ITS, nLSU, and EF-1a) with an alternative topology based on 
two mitochondrial loci (ATP6 and mtSSU) (Fig. 5). The nuclear 
phylogeny showed no obvious phylogeographic divergence 
within individual taxa. In contrast, the mitochondrial tree showed 
molecular divergence that unexpectedly followed patterns of 
geography rather than morphology. The case was particularly 
striking within the Australasian Rossbeevera species where 
we detected little differentiation between species based on the 
three nuclear loci. In contrast, the mitochondrial genes largely 


exhibited a pattern of isolation by distance with geographically 
distant specimens also being genetically distant. Similar pat- 
terns were apparent in several species of Japanese Rossbee- 
vera and Turmalinea (Fig. 5). Specimens of R. griseovelutina 
show considerable infraspecific divergence in the nuclear 
phylogeny. Specimens collected from western Honshu (the 
mainland of Japan) formed a single, homogenous clade based 
on the mitochondrial dataset whereas specimens from eastern 
Honshu and Amami-oshima Island (c. 600 km from Honshu) 
were genetically unique and more closely related to each other. 
In T. persicina, the nuclear phylogeny did not exhibit any obvi- 
ous phylogeographic pattern, but the mitochondrial phylogeny 
suggested that specimens from Kyushu were slightly divergent 
from other Honshu specimens. These findings indicate that the 
infraspecific lineages in the genus Rossbeevera (including the 
Australasian species complex) may not be reproductively iso- 
lated, and that the infraspecific divergences in the mitochondrial 
genes are correlated to phylogeographic divergence rather than 
morphological and ecological differences. 


Disagreement between the phylogenetic signals in nuclear vs 
mitochondrial loci have been documented for epigeous fungal 
genera and families such as the saprobic genus Sparassis 
(Polyporales) (Dai et al. 2006), ectomycorrhizal genus Strobilo- 
myces (Boletales) (Sato et al. 2007, Sato & Murakami 2008), 
and the Clavicipitaceae, which includes fungi with diverse 
trophic modes (Sung et al. 2007). However, studies of these 
fungi did not detect any obvious biogeographic pattern in the 
diversification of mitochondrial loci (ATP6 and/or COX3 genes). 
It is possible that the unique ecological constraints imposed on 
sequestrate fungi (e.g. belowground fruiting, lack of active spore 
dispersal) may be correlated with geographical divergence in 
the mitochondrial phylogeny, but more data and comparisons 
with additional sequestrate taxa and related epigeal taxa are 
clearly needed to resolve this question. 


The comparison between the nuclear and mitochondrial phylo- 
genies also revealed divergent topological patterns within 
infraspecific lineages of Rossbeevera. One R. griseovelutina 
specimen (TNS-F-36991) is resolved in the R. griseovelutina 
lineage 2 based on the nuclear markers (ITS, nLSU, EF-1a) 
but is grouped together with the members of the R. griseovelu- 
tina lineage 1 based on the mitochondrial markers (ATP6 and 
mtSSU) (Fig. 5). Since the other specimens of R. griseovelutina 
lineage 2 in the nuclear tree showed phylogenetic affinity in the 
mitochondrial tree, ILS between the two lineages is not likely. 
This conflict could be more reasonably explained by introgres- 
sion of mitochondrial genes because both R. griseovelutina 
lineage 1 and lineage 2 are found in western Honshu where 
lineage 1 occurs in evergreen Castanopsis forests and lineage 2 
is found in nearby deciduous Fagus-Quercus forests. It is pos- 
sible that distribution of the two lineages overlaps in some areas 
in western Honshu, and this could cause the introgression of 
mitochondrial genes. It is not obvious, however, why the two 
individuals of the lineages 1 and 2 collected from Amami-oshima 
Island and Eastern Honshu, which are geographically distant 
from one another (c. 1 200 km), have almost identical ATP6 
and mtSSU sequences. Given that the mitochondrial sequences 
of Rossbeevera tend to reflect geographic distance, transoce- 
anic spore dispersal as well as gene introgression would be 
necessary to explain these topologies. The precise pattern 
of migration cannot be concluded thus far due to insufficient 
samples. However, additional sampling throughout the region 
might help to clarify the process of transoceanic spore dispersal 
of truffles in the future. 


Several studies on the evolutionary biology of fungi and ani- 
mals have detected similar incongruence between nuclear 
and mitochondrial datasets. Dai et al. (2006) demonstrated 
that three Sparassis species distinct in the nuclear phylogeny 
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shared identical mitochondrial ATP6 sequences. The authors 
inferred that this was due to clonal inheritance within the partly 
overlapping populations. Similarly, Robinson et al. (2001) found 
remarkably low divergence in ATP6 of Agaricus bisporus iso- 
lates with different mitochondrial haplotypes but they did not 
mention the possibility of ILS or introgression. Peters et al. 
(2007) examined the mitochondrial paraphyly between gadwalls 
and falcated ducks and rejected the hypothesis of ILS based 
mainly on incongruence between the divergence time of the 
two species estimated by nuclear introns and mtDNA genes. 
Similar approaches and addition of sequences of R. griseo- 
velutina may offer further insights about evaluation of ILS and 
introgression in Rossbeevera. 


Although ILS seems unlikely in the case of the topological con- 
flict within R. griseovelutina, it seems probable that ILS does 
account for the topological incongruence within R. paracyanea 
(Fig. 5). Both infraspecific lineages of R. paracyanea were from 
the same region (western Honshu, Japan), and Fig. 5 shows 
the infraspecific divergence in the mitochondrial two-locus 
phylogeny is deeper compared to that in the nuclear three- 
locus phylogeny. Although we cannot conclude any historical 
event that affects the topological incongruence based only on 
the small number of R. paracyanea samples, the parsimoni- 
ous explanation is incongruence caused by ILS derived from 
mitochondrial polymorphysm within the ancestor of the extant 
R. paracyanea lineages. Considering these cases of putative 
ILS between the nuclear and mitochondrial loci and the likely 
ILS found among the nuclear single-locus trees, it seems that 
ILS is not a rare event in this group of fungi. Furthermore, the 
frequent polymorphism in the nuclear loci of the Australasian 
Rossbeevera ‘spp.’ may be another example of the ongoing 
process of lineage sorting (Fig. 1). 


Similar topological patterns to R. paracyanea were recovered 
in the relationship between R. cryptocyanea lineages 1 and 2 
(Fig. 5). In this case, however, both ILS and introgression are 
possible causes of the topological discordance because the 
pattern of divergence in the mitochondrial tree also reflected 
geographical distances among the specimens. The distance be- 
tween the two localities, Kakeroma and Amami-oshima Islands, 
is c. 15 km, whereas Ishigaki Island is c. 650 km distant from 
Amami-oshima Island, and the Kerama Gap, which is known as 
a faunal gap of many terrestrial animals, crosses between them 
(Ota 1998). Therefore, it is also possible that the two different 
infraspecific lineages of R. cryptocyanea occurred sympatri- 
cally and that this caused introgression of mitochondrial genes. 


Overall, the analysis of nuclear vs mitochondrial loci revealed 
that genetic distances in the mitochondrial tree generally cor- 
respond to geographic distance within each species-level line- 
age of Rossbeevera and Turmalinea. While some topological 
conflicts within Rossbeevera species may be the results of 
recent gene introgression, other conflicts strongly suggest the 
possibility of incomplete lineage sorting between closely related 
lineages. These examples indicate that careful, multifaceted 
characterization using both morphological and phylogeographi- 
cal methods with nuclear and mitochondrial datasets is the 
key to delimit boundaries between infrageneric taxa. It is also 
possible that these types of complex evolutionary histories 
may be frequently overlooked in phylogenetic studies of many 
other fungal taxa due to insufficient sampling of individuals and 
due to the lack of precise topological comparisons between 
phylogenies based on different loci. 


Minisatellite-like insertion in ITS2 as a DNA barcode marker 


The ITS region has recently been adopted as the official mo- 
lecular barcode for fungi (Chase & Fay 2009, Schoch et al. 
2012). Although most fungi in Basidiomycota have an ITS region 
(ITS1-5.8S-ITS2) that is approximately 500—650 nucleotides, 


PCR amplicons of ITS in Rossbeevera and Turmalinea species 
are remarkably large (e.g., ITS in T. persicina is about 1 100 
base pairs). The extra length of this region is due to the pres- 
ence of a minisatellite-like insertion within ITS2 (Table 2). ITS 
minisatellites are rare in fungi, but Den Bakker et al. (2004) 
reported that species of Leccinum, a species-rich bolete genus 
in the leccinoid clade, usually contain minisatellite-like insertions 
in ITS1. They demonstrated that the minisatellite-like insertions 
were highly divergent among infraspecific lineages of Leccinum 
and therefore not usable as a phylogenetic marker at or above 
species level. Orihara et al. (2012b) briefly mentioned the pres- 
ence of insertions in the ITS2 region of species in Octaviania 
subg. Fulvoglobus, but did not characterize the DNA sequences 
of the insertions. 


In this study we used the Automatic Barcode Gap Discovery 
(ABGD) software (Puillandre et al. 2012), to infer species boun- 
daries within Rossbeevera and Turmalinea. The barcode gap 
analysis and ML phylogeny clearly showed that the minisatellite- 
like insertions of Rossbeevera and Turmalinea were highly 
conserved at the subspecies-, species- and genus levels and 
that these minisatellites have higher resolution in discriminating 
the boundary between intergeneric-level, infraspecific-level, 
and interspecific-level divergence than the rest of the ITS 
region (Fig. 6). The Australasian Rossbeevera spp., however, 
were not partitioned at the species-level boundary estimated 
by the analysis (Fig. 6c). This suggests that these previously 
described species are likely geographic variants that may still 
be in the process of speciation, despite the fact that several 
specimens from different regions of Australasia show remark- 
able morphological divergence (Lebel et al. 2012a). 


Although the exact origin of these ITS insertions are unclear, 
it is likely that they originated within the leccinoid clade since 
no other group in the Boletales appear to have them. Although 
we have not characterized the structure and function of the 
minisatellite-like insertions, the insertion sequences of Ross- 
beevera and Turmalinea are highly informative and can thus 
be used efficiently for DNA barcoding of fruitbodies or environ- 
mental sequences. We caution, however, that the extended ITS 
length in Rossbeevera and Turmalinea may actually reduce the 
likelihood that these DNA sequences will be recovered from 
environmental samples since longer fragments are sometimes 
discriminated against in mixed pools of DNA (Huber et al. 2009, 
Bellemain et al. 2010). 


Conclusions 


Our finding of a new sequestrate genus, Turmalinea, reinforces 
the fact that the leccinoid clade is remarkably rich in sequestrate 
fungi. The nuclear three-locus phylogenies strongly supported 
the distinctness of the four new species and one new subspe- 
cies of Turmalinea as well as the two new species of Ross- 
beevera despite the fact that precise infrageneric relationships 
within Rossbeevera remain unclear. Topological comparison 
between nuclear and mitochondrial phylogenies documented 
significant infraspecific biogeographical diversification of two 
mitochondrial loci in the Australasian Rossbeevera species 
but very low differentiation in nuclear loci. Pairwise comparison 
among nuclear and mitochondrial trees further suggested that 
mtDNA introgression and ILS have occurred within multiple 
inter- and infraspecific lineages of Rossbeevera. Finally, we 
used the recently developed barcode gap analysis, ABGD, to 
demonstrate that the minisatellite-like insertion found in the 
ITS2 of Rossbeevera and Turmalinea is highly informative for 
identifying infra- and interspecific and intergeneric divergence. 
Although the ABGD approach has only rarely been used to 
delimit species of fungi, our results suggest that this method 
holds great promise for future studies of molecular identification 
in fungi, particularly since the ITS region has been adopted as 
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the official fungal barcode region (Chase & Fay 2009, Schoch 
et al. 2012). 


In a few lineages of Rossbeevera and Turmalinea, we were 
not able to include multiple infraspecific OTUs in the phylo- 
genetic analyses because these sequestrate fungi are mostly 
hypogeous and they are rarely collected. Moreover, we know 
of no DNA sequences from R. mucosa or R. bispora because 
all of the available specimens were collected too long ago and 
some of them are fragmentary. We hope that in the near future 
comprehensive taxon sampling of both epigeous and seques- 
trate leccinoid taxa will reveal more precise phylogenetic and 
phylogeographic relationships and evolutionary processes. 
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